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This  thesis  develops  an  analytical  method  for  determining  an  optimal  level 
of  specialization  and  optimal  task  allocation  for  a  maintenance  manpower  force. 
The  method  assumes  that  maintenance  tasks  are  generated  by  a  system  of  identical 
machines  which  experience  random  malfunctions  and  require  periodic  service.  The 
impact  of  alternative  manpower  structures  on  system  performance  is  evaluated  us¬ 
ing  a  queuing  network  model.  Markov  decision  analysis  is  employed  to  determine  an 
optimal  assignment  of  manpower  resources  to  pending  tasks  as  the  network  status 
varies  over  time.  A  linear  programming  algorithm  is  derived  to  enable  simultaneous 
optimization  of  specific  assignment  decisions  and  the  overall  manpower  structure. 
The  optimization  method  is  developed  and  demonstrated  through  a  simple  exam¬ 
ple,  but  the  dimensionality  issues  associated  with  larger  system  models  are  also 
addressed.  The  method  is  specifically  applied  to  the  problem  of  maximizing  mil¬ 
itary  aircraft  sortie  generation  subject  to  a  constraint  on  maintenance  manpower 
expenditure. 
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This  thesis  develops  an  analytical  method  for  determining  an  optimal  level 
of  specialization  and  optimal  task  allocation  for  a  maintenance  manpower  force. 
The  method  assumes  that  mcdntenance  tasks  are  generated  by  a  system  of  identical 
machines  w^hich  experience  random  malfunctions  and  require  periodic  service.  The 
impact  of  alternative  manpower  structures  on  system  performance  is  evaluated  us¬ 
ing  a  queuing  network  model.  Markov  decision  analysis  is  employed  to  determine  an 
optimal  assignment  of  manpower  resources  to  pending  tasks  as  the  network  status 
varies  over  time.  A  linear  programming  algorithm  is  derived  to  enable  simultaneous 
optimization  of  specific  assignment  decisions  and  the  overall  manpower  structure. 
The  optimization  method  is  developed  and  demonstrated  through  a  simple  exam¬ 
ple,  but  the  dimensionality  issues  associated  with  larger  system  models  are  also 
addressed.  The  method  is  specifically  applied  to  the  problem  of  maximizing  mil¬ 
itary  2drcraft  sortie  generation  subject  to  a  constraint  on  maintenance  manpower 
expenditure. 
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Chapter  1 


Introduction 


1.1  Motivation 

The  United  States  Air  Force  operates  an  inventory  of  over  3.000  tactical 
fighter  aircraft  at  locations  throughout  North  America,  Europe,  and  the  Pacific. 
These  aircraft  are  extremely  complex  machines  which  require  maintenance  support 
from  a  variety  of  specialized  personnel.  Due  to  pressures  for  peacetime  monetary 
efficiency,  fighter  units  and  support  forces  typically  operate  at  fixed  installations 
which  can  harbor  large  concentrations  of  aircraft.  These  centralized  forces  are  be¬ 
coming  increasingly  vulnerable  to  qualitative  improvements  in  the  aerial  bombard¬ 
ment  and  surface-to-surface  missile  capabilities  of  potential  wartime  adversaries. 
Consequently,  Air  Force  leadership  has  endorsed  a  deployment  strategy  which  “calls 
for  decentralized,  small-unit  autonomy  with  the  mobility  and  flexibility  to  survive 
and  sustain  dispersed  combat  operations”  [38,  p.  2].  Planners  are  now  considering 
unit  deployments  involving  very  small  numbers  of  aircraft.  Unfortunately,  such 
decentralization  could  substantially  increase  maintenance  manpower  requirements. 
One  study  h2LS  shown  that  dispersing  aircraft  to  small  unit  bases  could  increase 
needed  manpower  by  two-thirds  just  to  ensure  the  avsulability  of  minimum  crew 
sizes  for  all  specialized  tasks  [4].  Even  more  manpower  would  be  needed  to  limit 
queuing  delays  to  a  level  where  current  standards  of  operational  effectiveness  could 
be  maintained. 

The  apparent  need  for  more  manpower  competes  with  the  constraints  im- 
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posed  by  a  decreasing  budget  and  shrinking  supply  of  enlistment-aged  citizens.  An 
alternative  approach  under  consideration  is  a  restructuring  of  maintenance  special¬ 
ties  so  that  required  tasks  can  be  accomplished  by  fewer  personnel  having  a  wider 
range  of  skills.  The  higher  costs  of  recruiting,  training,  and  retaining  maintenance 
“generahsts”  may  be  more  than  offset  by  the  efficient  manpower  utilization  that 
could  be  preserved  under  dispersed  operations.  This  possibility  motivates  an  impor¬ 
tant  question:  what  specialization  strategy  will  maximize  operational  effectiveness 
for  a  given  level  of  manpower  expenditure? 

1.2  Problem  Statement 

The  general  problem  of  optimal  specialization  for  maintenance  manpower 
could  apply  to  myriad  military  and  industrial  concerns.  Many  enterprises  depend 
on  the  steady  operation  or  availability  of  complex  machiner}'.  For  example,  a 
transportation  entity  (airline,  trucking  firm,  overnight  delivery  service,  etc.)  cannot 
function  competitively  without  a  well-maintained  fleet  of  vehicles.  A  manufacturing 
facility  may  rely  critically  on  sophisticated  production  equipment.  Timely  main¬ 
tenance  of  these  machines  requires  skills  which  are  expensive  to  develop,  and  the 
personnel  possessing  these  skills  are  valuable  resources.  Nevertheless,  the  division 
of  labor  among  maintenance  personnel  can  often  reflect  traditional  organizational 
boundaries  rather  than  efficient  allocation  with  respect  to  the  overall  objectives  of 
the  parent  enterprise.  This  thesis  presents  a  general  method  for  determining  an 
optimal  specialization  strategy,  that  is,  for  optimizing  the  number  of  maintenance 
specialties  and  the  allocation  of  tasks  to  specialties.  The  optimization  method  ex¬ 
plicitly  recognizes  the  economic  tradeoff  between  the  lower  per-person  costs  and 
the  lower  average  utilizations  resulting  from  specialization. 

As  indicated  in  Figure  1.1,  determination  of  an  optimal  specialization  strat- 
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Figure  1.1;  Components  of  the  Manpower  Problem 


egy  requires  suboptimization  at  two  levels.  First,  an  optimal  strategy  cannot  1  e 
identified  without  determining  an  optimal  manpower  structure,  which  specifies  ex¬ 
act  quantities  for  each  type  of  specialist  in  the  manpower  force.  Further,  to  find 
an  optimal  structure,  it  is  necessary  to  determine  an  optimal  assignment  policy. 
The  assignment  policy  defines  comprehensive  decision  rules  for  dispatching  man¬ 
power  to  specific  pending  tasks  as  the  enterprise  operates  over  time.  The  number 
of  possible  manpower  structures  and  assignment  policies  can  be  quite  large  for  a 
particular  application,  so  determining  an  optimal  specialization  strategy  by  exhaus¬ 
tive  enumeration  may  be  computationally  intractable.  Enumerative  evaluation  of 
all  manpower  structures  and  assignment  policies  can  be  avoided  through  the  ana¬ 
lytical  method  developed  in  this  thesis. 
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1.3  Assumptions 

The  approach  pursued  in  developing  the  optimization  method  assumes  that 
the  modeled  enterprise  has  the  general  characteristics  listed  below.  Relaxation  of 
some  of  these  assumptions  is  considered  as  the  model  development  progresses. 

1.  The  enterprise  employs  a  fixed  inventory  of  N  machines  with  identical  po¬ 
tential  maintenance  requirements.  Any  machine  not  requiring  maintenance  is 
immediately  utilized  in  operations. 

2.  Operations  are  conducted  continuously,  so  the  goal  of  the  enterprise  is  to  max¬ 
imize  steady-state  performance.  This  is  achieved  by  maximizing  the  expected 
number  of  machines  operating  at  a  random  point  in  time. 

3.  A  finite  set  of  M  tasks  describes  all  possible  types  of  machine  maintenance. 
Each  task  m  €  {1,2,...,  M}  is  characterized  by  a  required  number  of  person¬ 
nel  Cm  and  a  known  average  completion  rate  /x^ .  The  time  required  to  perform 
each  task  is  represented  by  an  exponentially  distributed  random  variable  [33, 
pp.  201-209]. 

4.  A  required  maintenance  task  can  be  generated  by  a  corresponding  type  of  ran¬ 
dom  malfunction,  or  by  completion  of  a  time-limited  operational  activity  (e.g., 
refueling  after  each  aircraft  sortie).  The  operating  time  between  random  mal¬ 
functions  of  any  type  m  is  represented  by  an  exponentially  distributed  random 
variable  with  rate  Am  for  each  machine.  Each  type  of  random  malfunction  is 
independent  from  all  other  types,  and  all  machines  operate  independently. 

5.  The  duration  of  each  operational  activity  is  represented  by  an  exponentially 
distributed  random  variable  with  rate  /to-  Maintenance  tasks  generated  due  to 
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random  malfunctions  can  be  deferred  until  completion  of  the  next  operational 
activity,  so  a  machine  can  reach  a  condition  where  several  types  of  tasks  are 
required.  In  this  case,  the  required  tasks  can  be  performed  sequentially  or  in 
parallel,  depending  on  a  specified  compatibility  of  the  task  tj’^pes. 

6.  Each  type  of  maintenance  task  must  be  allocated  to  a  particular  type  of  main¬ 

tenance  specialty.  The  number  of  specialties  in  a  particular  specialization 
strategy  can  range  from  one  (one  specialty  for  all  tasks)  to  M  (a  separate  spe¬ 
cialty  for  each  task).  Thus,  the  number  of  specialties  Y  to  be  considered  across 
all  possible  strategies  cannot  exceed  an  upper  bound  of  =  2^^  —  1. 

7.  Continuous  manning  of  any  maintenance  specialty  y  G  {l,2,...,y}  incurs 
a  known  cost  of  Cy  units  per  specialist.  The  enterprise  is  constrained  to  an 
upper  bound  of  C  units  for  total  maintenance  manpower  expenditure.  Thus, 
any  feasible  manpower  structure  must  satisfy  the  constraint  CyXy  <  C 
where  Xy  is  the  number  of  type  y  specialists  employed  under  the  structure. 

These  assumptions  provide  a  framework  for  analyzing  the  interacting  ef¬ 
fects  of  machine  characteristics  and  manpower  availability  on  overall  system  per¬ 
formance.  At  any  point  in  time,  an  individual  machine  may  be  operating,  or  it  may 
require  a  finite  number  of  maintenance  tasks.  The  state  of  the  overall  system  can 
be  described  by  listing  the  number  of  machines  in  each  possible  condition.  If  the 
number  of  machines  requiring  a  particular  maintenance  task  exceeds  the  allocated 
number  of  specialists,  the  excess  machines  will  queue  for  an  available  specialist 
and  potential  operating  time  will  be  lost.  Thus,  each  maintenance  condition  could 
correspond  to  a  station  in  a  “queuing  network,”  with  the  stations  connected  by 
the  possible  movement  of  a  machine  from  one  condition  to  another.  Since  oper¬ 
ating  times  and  maintenance  times  are  modeled  with  exponential  distributions. 
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the  movement  of  machines  through  the  network  can  be  described  as  a  continuous¬ 
time  Markov  process  [37,  pp.  39-45].  Several  analytical  methods  can  be  applied  to 
evaluate  a  system  conforming  to  this  general  model. 

1.4  Literature  Review 

The  published  literature  generally  associated  with  the  disciplines  of  in¬ 
dustrial  engineering,  management  science,  and  operations  research  makes  scarce 
reference  to  the  specific  issue  of  optimal  manpower  specialization.  However,  the 
general  topic  of  resource  utilization  in  a  queuing  environment  is  well  developed, 
and  many  theoretical  contributions  can  be  apphed  to  the  manpower  problem.  In 
addition,  several  government  studies  have  quantitatively  addressed  some  manpower 
issues  associated  with  complex  military  systems.  Many  types  of  machines  with  de¬ 
fense  applications  (particularly  aircraft)  tend  to  have  a  high  ratio  of  maintenance 
man-hours  to  operating  hours,  and  the  resources  required  to  perform  maintenance 
are  expensive. 

1.4.1  Queuing  Network  Analysis 

The  high  requirements  and  costs  of  aircraft  maintenance  motivated  the  first 
published  investigations  of  queuing  network  behavior.  In  1954,  the  concept  of  a 
cyclic  queue  was  introduced  by  J.  Taylor  and  R.  R.  P.  Jackson  as  a  model  for 
relating  maintenance  rates  and  spare  engine  supply  to  the  availability  of  a  fleet  of 
aircraft  [36].  Taylor  and  Jackson  considered  the  movement  of  ?  finite  supply  of 
customers  (aircraft  engines)  through  a  series  of  service  stations  with  infinite  queue 
capacities  and  exponentially  distributed  service  times.  In  the  same  journal,  R.  R. 
P.  Jackson  published  another  article  which  addressed  a  similar  problem  with  an 
infinite  arrival  population  [19].  Jackson  theorized  that  the  output  distribution  of  a 


7 


service  station  was  identical  to  the  input  distribution.  This  result,  proved  by  Burke 
in  1956  [7],  permitted  independent  treatment  of  each  queue  in  a  series.  In  1957, 
J.  R.  Jackson  published  Networks  of  Waiting  Lines,  which  formalized  the  earlier 
results  of  R.  R.  P.  Jackson  and  extended  them  to  include  systems  of  queues  in 
which  transitions  of  customers  from  one  queue  to  another  at  service  completion  are 
random  events  [20].  Such  systems  are  now  often  referred  to  as  “Jackson  networks.” 

In  the  original  Jackson  network  model,  customers  arrive  from  an  infi¬ 
nite  external  source  according  to  a  Poisson  process.  Each  service  station  con¬ 
tains  one  or  more  parallel  servers,  and  all  service  times  are  exponentially  dis¬ 
tributed.  Any  state  of  a  system  with  Z  +  1  stations  can  be  represented  by  a 
vector  n  =  (7io,n], . . .  ,nz)  where  >  0  is  the  number  of  customers  at  each  sta¬ 
tion  z  €  {0,1,...,Z}.  The  equilibrium  probability  that  the  network  is  in  state  n 
is  denoted  by  P{n)  =  P(no,ni, . .  .,nz),  and  the  marginal  probability  of  finding  n^ 
customers  at  a  station  z  is  denoted  by  ^2(71^).  Jackson  proved  that  the  equilib¬ 
rium  probability  of  a  given  state  n  can  be  factored  into  the  product  of  each  of  the 
marginal  distributions;  that  is, 

z 

P{no,ni,...,nz)  =  Yipzin,)  (1.1) 

z=0 

This  “product  form”  relationship,  coupled  with  Burke’s  earlier  result,  greatly  sim¬ 
plified  the  analysis  of  many  queuing  systems. 

In  1967,  Gordon  and  Newell  extended  the  product  form  resiJt  to  obtain  the 
equilibrium  state  probabilities  for  closed  Markovian  networks  [14].  They  studied  a 
system  with  a  finite  number  of  customers  N  cycling  through  a  network  of  stations 
with  multiple  parallel  servers.  The  state  of  the  system  is  again  described  by  a 
vector  n  =  {no,ni, . . .  ,nz)i  but  now  an  additional  constraint  is  imposed;  that  is. 


8 


IZf=o  Gordon  and  Newell  proved  that  the  equilibrium  probability  of  state  n 

also  satisfies  a  product  form 

1  ^ 

P(no,ni,...,n^)  =  JJ/,(n,)  (1.2) 

where  each  function  depends  only  on  the  characteristics  of  the  zth  station  and 
T{N)  is  a  normalization  constajit  chosen  to  make  all  the  feasible  state  probabil¬ 
ities  sum  to  one.  Although  the  product  form  developed  by  Gordon  and  Newell 
was  easily  expressed,  direct  calculation  of  the  state  probabilities  was  often  com¬ 
putationally  expensive  and  inaccurate  because  of  the  large  number  of  operations 
required  to  evaluate  the  normalization  constant.  In  1973,  J.  P.  Buzen  developed 
an  efficient  algorithm  for  evaluating  r(iV)  [8].  Buzen  also  studied  methods  for  ob¬ 
taining  the  performance  measures  of  a  queuing  system  as  a  simple  function  of  the 
model  parameters  and  the  normalization  constant. 

Certain  types  of  queuing  networks  do  not  conform  to  the  product  form 
model.  In  1978,  P.  J.  Denning  surveyed  an  “operational  approach”  to  queuing  sys¬ 
tem  analysis  and  thoroughly  studied  the  conditions  required  for  solution  of  queuing 
network  models.  The  first  condition  needed  for  determination  of  equilibrium  state 
probabilities  is  “flow  balancing,”  meaning  that  the  transition  rates  in  and  out  of 
each  system  state  must  be  balanced.  This  condition  permits  construction  of  a  set 
of  global  balance  equations  (one  equation  for  each  state)  which  can  be  solved  simul¬ 
taneously  by  traditional  methods.  Unfortunately,  even  modest-sized  networks  can 
generate  a  large  number  of  simultaneous  equations.  The  product  form  approach 
circumvents  this  problem,  but  requires  that  two  additional  conditions  be  met: 

1.  Multiple  customer  arrivals  and  departures  are  not  observed,  so  that  the  rate  at 
which  the  system  enters  or  leaves  a  state  depends  oidy  on  the  rate  of  customer 
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flow  between  stations  (“one-step  behavior”). 

2.  The  flow  rate  out  of  a  station  depends  only  on  the  station’s  queue  length,  and 
not  on  how  customers  are  distributed  elsewhere  in  the  system  (“homogene¬ 
ity”). 

Under  these  circumstances,  the  global  balance  equations  decompose  into  a  set  of 
local  balance  equations  and  the  product  form  result  can  be  applied.  However,  not 
all  queuing  systems  satisfy  the  required  conditions.  For  example,  the  homogeneity 
condition  is  violated  by  systems  which  involve  some  sort  of  “blocking,”  where  an 
event  at  one  station  can  prevent  another  station  from  serving  customers.  This  can 
occur  if  resource  conflicts  exist,  such  as  when  queues  must  share  servers.  Unfor¬ 
tunately,  this  is  precisely  the  situation  that  arises  when  several  maintenance  tasks 
require  one  type  of  resource  and  total  task  requirements  exceed  the  number  of 
available  resources.  Thus,  the  manpower  specialization  problem  requires  solution 
techniques  which  do  not  rely  on  product  form  methods. 

Researchers  concede  a  lack  of  progress  in  obtaining  exact  solutions  for  large 
queuing  networks  that  do  not  have  product  form  state  probabilities  [29].  Two  prin¬ 
cipal  approximation  methods  for  attacking  these  networks  have  emerged:  diffusion 
and  decomposition.  The  diffusion  approximation  attempts  to  reduce  the  com¬ 
plexity  of  the  global  balance  equations  by  treating  individual  queue  lengths  as 
continuous  random  variables  [21].  This  technique  can  be  particularly  useful  for  han¬ 
dling  nonexponential  service  time  distributions  and  for  finding  transient  solutions 
[22].  Unfortunately,  diffusion  methods  can  be  difficult  to  implement  and  are  based 
on  an  assumption  of  “heavy  traffic”  which  would  always  favor  highly  specialized 
manpower  structures.  Thus,  approximation  methods  based  on  diffusion  are  less 
applicable  to  the  manpower  specialization  problem  than  approximation  methods 
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based  on  decomposition. 

The  decomposition  approach  was  introduced  in  1975  bj"  P.  J.  Courtois  as 
a  tool  for  performance  evaluation  of  complex  computer  systems  [11].  An  example 
of  this  approach  might  consist  of  analyzing  a  subsystem  in  isolation  and  then  re¬ 
placing  the  subsystem  with  a  single  composite  server  which  imitates  the  behavior 
of  the  originally  isolated  subsystem.  Courtois  reported  that  such  an  approximation 
is  accurate  if  the  rate  of  interaction  within  the  subsystem  is  substantially  higher 
than  the  rate  of  interaction  between  the  subsystem  and  the  rest  of  the  overall 
network.  The  decomposition  approach  was  further  developed  by  Chandy,  Herzog, 
and  Woo,  who  applied  Norton’s  theorem  from  electrical  circuit  theory  to  queuing 
networks  which  obey  local  balance  [9].  Chandy  et  al.  demonstrated  that,  given 
a  network  with  Z  -h  1  stations,  it  is  possible  to  replace  Z  of  the  stations  with  a 
single  composite  server  having  load-dependent  service  rates.  They  also  presented 
an  approximate  procedure  for  modeling  this  “complementary  queue”  in  a  network 
which  does  not  obey  local  balance.  These  techniques  permit  construction  of  ap¬ 
proximate  models  of  otherwise  intractable  systems  by  defining  stations  which  are 
“flow  equivalent”  to  major  portions  of  an  original  network.  This  notion  of  flow 
equivalence  has  significant  applicability  to  the  problem  of  determining  an  optimal 
manpower  specialization  strategy  for  a  large  system. 

1.4.2  Military  Manpower  Analysis 

The  high  maintenance  requirements  of  complex  military  systems  suggest 
that  maintenance  manpower  would  be  a  major  consideration  in  system  design. 
However,  quantitative  methods  have  not  always  been  applied  to  thoroughly  ad¬ 
dress  manpower  issues.  In  1960,  a  research  memorandum  produced  by  the  RAND 
Corporation  for  the  U.  S.  Air  Force  reported: 
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Manpower  planning,  until  recently,  has  played  a  rather  passive  role  in 
the  research  and  development  of  future  Air  Force  weapon  and  support 
systems.  Plans  for  the  best  utilization  of  people  in  new  systems  have 
been  developed  late  in  the  R&D  cycle  and  usually  after  decisions  have 
been  made  about  hardware  and  basic  operational  characteristics.  [16, 

p.  1] 

The  RAND  memorandum  described  a  general  approach  for  analyzing  manpower  re¬ 
quirements,  but  did  not  provide  specific  quantitative  techniques.  In  the  years  since 
this  early  report,  manpower  issues  have  been  studied  extensively,  but  almost  exclu¬ 
sively  through  models  which  employ  discrete  event  simulation  of  logistic  systems 
[30].  One  notable  exception  is  the  DYNA-METRIC  model,  which  is  au  analytical 
tool  used  primarily  to  measure  the  impact  of  spare  parts  management  on  aircraft 
sortie  generation  [17,  34].  In  current  manpower  analysis,  the  tool  most  frequently 
used  is  a  widely  accepted  computer  simulation  known  as  the  Logistics  Composite 
Model  (LCOM). 

Simulation  models  can  be  constructed  to  capture  an  arbitrary  level  of  detail, 
but  the  price  for  this  advantage  includes  uncertain  accuracy,  substantial  expense, 
and  a  lack  of  analytical  insight  into  causal  relationships.  Simulation  is  particularly 
poorly  suited  for  optimization  problems,  since  the  random  nature  of  results  hinders 
iterative  convergence  on  an  optimal  combination  of  many  input  variables.  Math¬ 
ematical  queuing  models  have  occasionally  been  employed  in  manpower  studies, 
but  usually  in  analyzing  particular  system  components  rather  than  the  interactions 
between  all  maintenance  requirements  and  available  manpower.  The  manpower 
specialization  problem  addressed  in  this  thesis  involves  a  broader  context  than  has 
thus  far  been  addressed  with  mathematical  methods. 

In  1964,  H.  R.  Barton  and  others  developed  an  early  analytical  method  for 
investigating  the  logistical  requirements  of  alternative  system  designs  [1,  32].  The 
method  employed  queuing  tables  to  analyze  tradeoffs  between  manpower  levels. 
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spare  parts  inventory,  and  machine  downtime.  Machine  components  were  consid¬ 
ered  separately  for  specified  levels  of  service  demand  that  were  externally  deter¬ 
mined.  This  general  approach  typified  subsequent  applications  of  queuing  theory 
to  military  logistics  problems,  including  the  DYNA-METRIC  model  and  its  prede¬ 
cessors  [17,  p.  11]. 

In  1975,  J.  R.  Phelan  employed  simple  queuing  network  analysis  to  predict 
the  impact  of  maintenance  manpower  levels  on  aircraft  operational  effectiveness. 
Different  types  of  maintenance  requirements  were  modeled  in  separate  cycles,  and 
interactions  were  explicitly  ignored  [31,  p.  23].  Interestingly,  the  computational 
methods  used  in  this  analysis  relied  on  global  balance,  even  though  product  form 
solution  techniques  could  have  been  applied.  The  effort  also  presented  some  jus¬ 
tification  for  considering  manpower  issues  separately  from  other  dimensions  of  the 
overall  logistics  system  (spare  parts,  etc.). 

Beginning  in  the  mid  1980s,  the  specific  problem  of  optimal  task  allocation 
for  aircraft  maintenance  specialties  became  the  subject  of  keen  interest.  This  inter¬ 
est  was  partially  motivated  by  increased  equipment  sophistication  and  maintenance 
complexity,  which  raised  the  aptitudes  and  skills  required  of  maintenance  person¬ 
nel.  Competing  demands  for  the  same  aptitudes  and  skills  in  other  occupations 
(military  and  industrial),  coupled  with  projected  demographic  changes,  presented 
serious  manpower  challenges.  These  challenges  were  amplified  by  an  institutional 
goal  to  improve  aircraft  survivability  through  dispersal,  an  operational  concept  that 
would  stretch  maintenance  resources  even  further. 

In  1983,  M.  Berman  and  C.  Batten  used  the  TSAR  simulation  model  (The¬ 
ater  Simulation  of  Airbase  Resources)  to  estimate  the  number  of  aircraft  sorties 
that  could  be  generated  under  differing  degrees  of  dispersal  [4].  Variables  consid¬ 
ered  in  the  analysis  included  manpower,  spare  parts,  and  aircraft  reliability  and 
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maintainability  characteristics.  While  Berman  and  Batten  did  not  explicitly  con¬ 
sider  maintenance  costs,  they  concluded  that  specialty  consolidation  could  improve 
operational  performance  or  reduce  manpower  requirements. 

Simulation  analysis  was  also  employed  by  C.  H.  Shipman  in  1985  to  evalu¬ 
ate  the  potential  manpower  savings  realizable  through  specialty  consolidation.  He 
employed  the  Logistics  Composite  Model  to  simulate  ground  attack  fighter  opera¬ 
tions.  Shipman  found  that  a  dispersed  unit  of  eighteen  aircraft  could  be  maintained 
by  27%  fewer  technicians  through  “minor”  specialty  consolidation  (combining  two 
or  three  specialties  into  one),  and  by  37%  fewer  technicians  through  “major”  spe¬ 
cialty  consolidation  (combining  six  specialties  into  one).  Shipman  recognized  the 
existence  of  “tradeoffs  between  the  manpower  savings  of  combining  specialties  and 
the  additional  skill/ training  requirements  these  combinations  generate”  [35,  p.  16]. 
He  recommended  that  these  tradeoffs  be  studied  to  determine  an  optimal  level  of 
consolidation. 

In  1986,  G.  A.  Gotz  and  R.  E.  Stanton  considered  the  impact  of  main¬ 
tenance  specialty  cross-training  on  operational  effectiveness  [15].  They  simulated 
a  simple  aircraft  operation  that  involved  two  types  of  repairable  components  and 
two  maintenance  specialties.  Gotz  and  Stanton  concluded  that  cross-training  is 
particularly  important  if  wartime  failure  rates  and  repair  rates  differ  significantly 
from  pre-war  expectations.  They  also  pointed  out  that,  when  maintenance  special¬ 
ists  with  multiple  skills  are  introduced,  decision  rules  must  also  be  implemented  to 
specify  which  malfunctioning  components  will  be  repaired  first. 

The  Air  Force  recently  engaged  in  a  major  initiative  to  broadly  address 
maintenance  manpower  issues.  The  initiative,  called  “Rivet  Workforce,”  succeeded 
in  generating  some  institutional  momentum  toward  specialty  consobdation.  The  re¬ 
sulting  need  for  an  analytical  framework  to  evaluate  particular  consolidation  strate- 
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gies  motivated  a  recent  study  effort  called  SUMMA,  or  Small  Unit  Maintenance 
Manpower  Analysis  !24.  28,  38].  A  major  product  of  the  SUMMA  project  is  a 
computerized  decision  aid  to  assist  in  the  derivation  of  optimal  task  allocations. 
The  computer  model  employs  a  practical  analytical  method  w'hich  is  based  on  as¬ 
sumptions  that  all  sortie  durations  and  maintenance  times  are  deterministic,  all 
aircraft  fly  sorties  and  receive  maintenance  in  “batches,”  and  all  aircraft  mainte¬ 
nance  tasks  are  performed  in  series.  The  task  allocations  and  manpower  estimates 
derived  through  the  analytical  method  are  refined  through  iterative  simulation  ex¬ 
periments.  While  the  complexity  of  the  overall  maintenance  system  guarantees  a 
significant  role  for  simulation  modeling,  substantial  insight  can  be  gained  through 
additional  mathematical  treatment  of  the  problem.  The  system  structure  suggests 
a  meaningful  role  for  queuing  network  analysis. 


Chapter  2 


Model  Development 

Development  of  a  manpower  specialization  model  using  queuing  network 
analysis  is  perhaps  best  illustrated  through  a  simple  example.  In  this  chapter, 
global  balance  is  used  to  obtain  an  optim2d  specialization  strategy  for  a  smaU 
commercial  enterprise  called  “Mike’s  Flying  Club.” 

2.1  Problem  Formulation 

Mike  operates  an  around-the-clock  flying  service  at  a  municipal  airport.  His 
rented  facility  will  support  only  two  aircraft  {N  =  2),  and  demand  for  flying  sorties 
is  high  enough  that  both  machines  can  be  employed  whenever  they  are  available. 
Mike  attributes  this  high  demand  to  his  reasonable  fees,  so  he  is  very  interested  in 
controlling  his  expenses.  He  is  particularly  concerned  about  the  high  costs  of  good 
aircraft  mechanics  and  is  determined  to  keep  his  maintenance  payroll  below  $100 
per  hour  (C  =  100). 

In  the  course  of  day-to-day  operations,  Mike’s  airplanes  can  require  any  of 
three  different  types  of  maintenance  tasks: 

1.  Routine  “turn-around”  maintenance,  such  as  refueling  and  routine  inspection, 
which  is  always  required  between  sorties.  This  task  is  completed  at  an  hourly 
rate  fix  =  1.0  by  a  single  mechanic  (aj  =  1). 

2.  Unplanned  “airframe”  maintenance,  which  is  required  when  an  aircraft  returns 
from  a  sortie  with  a  reported  malfunction  not  related  to  the  aircraft  engine. 
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This  task  is  accomplished  at  rate  fi2  =  -25,  and  also  requires  one  mechanic 

(02  =  1). 

3.  Unplanned  “engine”  maintenance,  which  is  required  when  an  aircraft  returns 
with  reported  engine  problems.  The  service  rate  for  this  task  is  /i3  =  .5,  and 
two  mechanics  are  required  (03  =  2). 

Occasionally,  a  returning  aircraft  wUl  require  both  airframe  and  engine  mainte¬ 
nance.  In  this  case,  the  two  types  of  repair  can  be  initiated  simultaneously,  provided 
sufficient  mechanics  are  available.  Turn-around  maintenance  is  never  initiated  until 
an  aircraft  completes  all  unplanned  airframe  and  engine  maintenance.  These  rela¬ 
tionships  are  illustrated  in  the  network  representation  shown  in  Figure  2.1.  Each 
station  in  the  network  can  be  characterized  by  a  set  of  pending  tasks  and  a  set 
of  eligible  tasks  E^.  Pending  tasks  are  those  which  must  be  accomplished  before  a 
machine  at  the  station  will  become  operational.  Eligible  tasks  are  those  which  can 
be  initiated  immediately  if  qualified  manpower  is  available. 

The  number  of  states  realizable  by  the  complete  system  is  determined  by 
the  number  of  ways  in  which  N  =  2  aircraft  can  occupy  Z  -I-  1  =  5  stations.  This 
is  a  standard  occupancy  problem  [13,  p.  36],  so  the  total  number  of  states  I  can 
be  calculated  as  =  (4)  =  15.  The  system  will  transition  to  a  new  state 

whenever  an  aircraft  completes  an  activity  (sortie  or  maintenance  task)  and  moves 
to  a  new  station.  Since  all  sortie  durations  and  maintenance  times  are  modeled  with 
exponential  distributions,  the  system  is  Markovian  and  exhibits  one-step  behavior. 
Each  feasible  transition  from  one  state  i  to  another  state  j  can  be  characterized  by  a 
losing  station  z~{ij)  and  a  gaining  station  z'*'{ij).  This  information  is  portrayed  in 
Table  2.1,  along  with  an  index  7(17)  which  identifies  the  activity  which  an  aircraft 
completes  at  transition.  A  value  of  =  0  implies  completion  of  an  aircraft 
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0:  aircraft  operating 

1:  turn-around  maintenance 

/>.  =  {!},£,  =  {1} 

2:  airframe  maintenance 

P2  =  {1,2},E2  =  {2} 

3;  engine  maintenance 
P,  =  {1,2},E,  =  {Z} 

4:  airframe  and  engine 
maintenance 

P,  =  {l,2,3},f;,  =  {2,3} 


Figure  2.1:  Possible  Conditions  for  Individual  Aircraft 


sortie,  while  a  value  of  'y(ij)  €  {1,2,3}  implies  completion  of  a  corresponding  type 
of  maintenance  task. 

The  transition  rates  between  system  states  will  depend  on  the  sortie  rate, 
maintenance  rates,  and  the  number  of  personnel  available  to  perform  the  various 
types  of  maintenance.  It  is  also  necessary  to  translate  aircraft  malfunction  rates 
into  network  routing  probabilities  for  aircraft  as  they  complete  sorties.  Aircraft 
maintenance  records  for  Mike’s  operation  reveal  average  airframe  and  engine  mal¬ 
function  rates  of  Aj  =  .2  and  A3  =  .25  respectively.  Flight  log-books  indicate  that 
an  average  sortie  has  a  duration  of  2.0  hours,  yielding  an  average  sortie  completion 
rate  of  fio  =  .5.  For  any  aircraft  beginning  a  sortie,  let  To  be  the  sortie  completion 
time,  let  T2  be  the  time  until  the  next  air&ame  malfunction,  and  let  T3  be  the  time 
until  the  next  engine  malfunction.  Routing  probabilities  92  for  z  E  {1,2, 3, 4}  can 
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be  computed  as  follows: 


9i 


P{To  <T2,To<  T3) 

P{Tu  <  min(T2,r3)} 

fio _ -5 

fiQ  +  A2  +  A3  .5  -h  .25  +  .2 


=  .5263 


92  =  P{To  >  T2,To  <  T3) 

=  P{To  <  T3)  -  P{To  <  T2,To  <  T3) 

flQ  .5 


/io  +  A^ 


-  9i 


.5  +  .25 


-  .5263  =  .1404 


93  =  P{Tq  <  T2.,To  >  T3) 

=  P{To<T2)-P{To<T2,To<T3) 

flQ  .5 


fi{]  +  A' 


-91  = 


.5  +  .2 


-  .5263  =  .1880 


94  =  P{To  >  TiiTq  >  ^3) 

=  1  -  P(To  <  T2,  To  <  T3)  -  P{To  >  T2,  To  <  T3)  -  P{To  <  T2,To  >  T3) 
=  1  -  -  92  -  93  =  1  -  .5263  -  .1404  -  .1880  =  .1454 


These  computations  can  be  generalized  for  an  arbitrary  network  with  any  number 
of  stations.  First,  let  A  be  the  set  of  tasks  which  are  always  required  between 
operational  activities.  Clearly,  any  station  whose  pending  tasks  do  not  include  all 
tasks  in  A  can  not  be  entered  directly  from  station  0.  Thus,  =  0  V  z  :  j4  g  P^. 
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For  all  other  stations, 


<lz 


P{To<Tm'^m^P,\jA)-  Y.  9c 


_ /fo _ 

"t"  ^dtn^PzUA 


E  ?c 

C-PcCPz 


(2.1) 


The  routing  probabilities  are  well  suited  for  sequential  computation  if  the  main¬ 
tenance  stations  are  ordered  by  increasing  numbers  of  tasks  in  their  pending  task 
sets. 

Once  the  routing  probabilities  are  computed,  global  balance  equations  can 
be  constructed  to  compute  the  equilibrium  state  probabilities  tt,  resulting  from  a 
particular  manpower  assignment  pobcy.  To  facilitate  a  general  representation  of  an 
assignment  pobcy,  it  is  useful  to  define  a  variable  s,j  as  the  number  of  mechanics 
assigned  to  perform  each  transition  task  7(1^)  6  {1,2,...,M}.  The  maximum 
value  of  each  s,j  is  Hmited  by  the  number  of  machines  requiring  maintenance; 
i.e.,  Sij  <  V  i,j  €  {1,2,...,/}.  The  values  are  also  constrained  by 

manpower  availabibty;  i’C;  Lj:-r(ij)eQy  <  ®v  V  t  G  {1,2, . . .  ,/},i/  e  {l,2,...,y}. 
Since  the  system  is  Markovian,  the  rate  of  a  transition  task  completion  can  clearly 
be  represented  as  It  is  also  clear  that  this  expression  must  reduce  to 

M’y(ij)  when  task  'y(ij)  is  the  only  one  pending  in  the  system,  since  there  would  be  no 
potential  for  a  resource  conflict.  Using  these  facts  and  the  information  presented  in 
Table  2.1,  global  balance  equations  for  each  state  of  the  example  can  be  constructed 
as  follows  (rate  out=rate  in): 


7ri2^0  =  TTz/il 


(2.2) 


+  Mo)  —  ’ri2gi^  -f  ^3^2  -f  ^4^3  -f  7r6562Mi 


(2.3) 


21 


7r3(M2  +  f^o)  —  ^i292Mo  +  ’*’5[“^]/^3  +  Tr7S73fl\ 


’T-lC/ia  +  A*o)  —  ’'■]2g3/iO  +  7f5554/i2  +  ’^r^SA^'l 


(2.4) 

(2.5) 


'’'’5{[^]/^3  +  ^54/^2  +  Mo}  —  7ri2g4Mo  +  TTgSgsMl 
X6562M1  =  ’'■291M0  -f  7r7S76M2  +  ’'‘8[-^]M3 


(2.6) 

(2.7) 


’'■7{573Mi  +  ^76M2}  =  ’r292Mo  +  +  ’'’gf— ]m3  +  ’rio5i0  7M2  +  ’'■ll[~y~]M3  (2.8) 


7r8{'S83Ml  +  [“^]/^3}  =  ’r293Mu  +  ’T^^iMo  +  7rgag8M2  +  TTiiSn  8M2  +  ’’■i3[~;^]M3  (2.9) 


^138i 


^9{^95Mi  +  [“^]/^3  +  ^98M2}  —  ^294MO  +  ^sQlf^O  +  7ri2"Sl2  9M2  +  ’’’mI  2  "  (2.10) 


r'SH9i 


I  r^l210, 

’'■io5i0  7M2  —  ’'■392M0  +  ’’’ni — 2 


(2.11) 


’ru{[-^]M3  +  -sii  8M2}  =  ’'■393M0  +  ’f492Mo  +  7712812 11M2  +  — ]/^3  (2.12) 


^14  11  • 


7ri2{5i2  9M2  +  [~^~]M3  +  ^12  nM2}  =  7r3g4Mo  +  7r5g2Mo  +  7ri5[  -  ]m3  (2.13) 


^15  ]2i 


r^l38i 

7ri3[— ^JM3  —  '^AQsfJ’O  +  7714814  13M2 


(2.14) 


’'’14{[^~]M3  +  [~^~]M3  +  •*14  13M2}  =  774g4Mo  +  +  77i58i5  14M2  (2.15) 


r*14  11  • 


’‘■l5{[  +  81514M2}  —  775q'4Mo 


(2.16) 


Since  flow  must  be  balanced  for  each  state,  any  one  of  the  above  equations  will 
be  redundant.  A  solution  for  the  state  probabilities  can  be  realized  by  arbitrarily 
deleting  Equation  2.16  and  replacing  it  with  the  normalizing  constraint 


(2.17) 


Once  the  TTj  are  obtained  by  simultaneous  solution  of  the  balance  equations,  the 
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expected  number  of  aircraft  operating  at  a  random  point  in  time  can  be  calculated 
as 


E{no)  =  ^  Tio^x,  =  2xi  +  7r2  +  TTa  +  TT^  +  TTs  (2.18) 

«=i 

Mike’s  expected  sortie  generation  rate  (sorties  per  aircraft  per  day)  can  be  com¬ 
puted  as 


R  —  24/^0 


N 


=  6£(no) 


(2.19) 


The  sortie  generation  rate  is  a  standard  measure  of  effectiveness  for  aircraft  oper¬ 
ations. 

It  should  be  emphasized  that  the  balance  equations  require  specific  input 
values  Sij  for  the  number  of  mechanics  designated  to  perform  specific  tasks  at 
specific  stations  in  each  system  state.  Twenty-two  of  these  values  appear  in  Equa¬ 
tions  2.2  through  2.16,  and  each  must  be  treated  as  a  variable  in  optimizing  system 
performance.  This  is  true  because  system  performance  will  be  determined  not  only 
by  a  prespecified  manpower  structure,  but  also  by  Mike’s  decisions  on  how  the 
manpower  will  be  utilized  when  each  system  state  is  encountered.  For  example, 
suppose  Mike  currently  employs  three  types  of  mechanics,  specializing  in  each  of  the 
three  different  types  of  maintenance  tasks  (maximum  specialization).  Specifically, 
suppose  that  Mike’s  payroll  includes  two  turn-around  mechanics  (zi  =  2),  one  air¬ 
frame  mechanic  (z2  =  !)>  and  two  engine  mechanics  (za  =  2).  When  the  system  is 
in  state  12,  so  that  n  =  (0, 0, 1,0, 1),  the  single  airframe  mechanic  could  be  assigned 
to  perform  an  airframe  task  on  the  aircraft  in  condition  2  (si2  9  =  l,5i2ii  =  0)  or 
the  aircraft  in  condition  4  (si29  =  0,S]2ii  =  1).  A  similar  conflict  exists  for  engine 
mechanics  in  state  14,  resulting  in  a  total  of  2  x  2  =  4  assignment  policies.  The 
potential  for  this  kind  of  resource  conflict  increases  if  specialization  is  reduced  so 
that  a  mechanic  can  perform  more  than  one  type  of  task.  Thus,  an  optimization 
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method  must  consider  not  only  an  overall  manpower  structure,  but  also  a  com¬ 
prehensive  policy  for  assigning  manpower  to  specific  tasks  when  resource  conflicts 
exist. 


2.2  Optimal  Specialization  Strategy 

Currently,  Mike  employs  shifts  of  turn-around,  airframe,  and  engine  me¬ 
chanics  at  hourly  wages  of  $10,  $20,  and  $25,  respectively.  He  has  observed  that 
mechanics  are  sometimes  idle  because  their  special  skills  don’t  always  correspond 
with  existing  maintenance  requirements.  Mike  knows  he  could  improve  manpower 
utilization  by  hiring  mechanics  who  are  qualified  in  multiple  types  of  maintenance, 
but  he  would  have  to  pay  significantly  higher  wages.  The  hourly  wage  for  each 
available  type  of  mechanic  y  €  {1,2, ...,5}  is  shown  in  Table  2.2.  Also  shown  is 
the  set  of  tasks  Qy  for  which  each  type  of  mechanic  would  be  qualified. 


Table  2.2:  Costs  and  Skills  of  Available  Manpower 


mm 

2 

3 

4 

5 

^9 

$10 

$20 

$25 

$30 

$33 

{1} 

Esm 

IBl 

Mike  wants  to  determine  the  maintenance  manpower  structure  (and  hence, 
the  specialization  strategy)  that  will  maximize  his  sortie  generation  rate  without 
exceeding  his  constrained  level  of  manpower  expenditure.  A  particular  manpower 
structure  can  be  represented  by  a  vector  x  =  (xj,  *2?  ••  •■>  *5)?  where  each  Xy  is  the 
number  of  type  y  mechanics  employed.  An  acceptable  structure  will  satisfy  the 
following  conditions: 
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1.  Total  hourly  manpower  expenditure  will  not  exceed  the  fixed  budget;  i.e., 
Ej=i  CyXy  <  100. 

2.  Each  task  type  will  be  assigned  to  only  one  type  of  mechanic;  i.e., 
if  Qy  n  Qy^  /  0,  then  XyXyi  =  0  V  e  {1,2, . . .  ,5},  y  ^  y' . 

3.  Sufficient  mechanics  will  be  qualified  to  perform  each  type  of  task;  i.e., 

Xy  >  TCliS.m^Qy{am)  V  j/  €  {1, 2, .  .  .  ,  5}. 

4.  Mechanics  of  each  type  will  have  a  potential  for  simultaneous  utilization;  i.e., 
Xy  <  maX;,g{i.2,3,4}(Em€£*nQ„  a„iV)  V  2/  €  {1, 2, . . . ,  5}. 

Table  2.3  lists  all  manpower  structures  which  satisfy  the  above  constraints  and  re¬ 
quire  a  total  expenditure  as  close  to  $100  as  possible  without  exceeding  this  limit. 
Let  (?  =  5  be  the  total  number  of  candidate  manpower  structures  and  =  3  be 
the  total  number  of  specialization  strategies.  For  each  structure  g  €  {1,2,...,Ct} 
and  corresponding  strategy  h  €  {1,2, . . . ,  if},  the  table  lists  the  number  of  poli¬ 
cies  generated  by  all  combinations  of  manpower  assignment  decisions.  To  find 
an  optimal  structure  by  exhaustive  enumeration,  it  wotild  be  necessary  to  solve 
4-l-2-l-24-fl2-f 36=78  sets  of  1  simultaneous  equations. 

Fortunately,  enumeration  can  be  avoided  by  employing  a  continuous-time 
Markov  decision  model  [18,  pp.  92-114].  Development  of  the  model  begins  with  the 
generation  of  a  list  of  manpower  assignment  decisions  which  are  possible  for  each 
state  across  an  unrestricted  range  of  manpower  structures.  Table  2.4  displays  all 
of  the  possible  decisions  k  for  the  “Mike’s  Flying  Club”  example.  For  each  specific 
structure  g  and  system  state  t,  a  subset  Dgi  of  nondominated  feasible  decisions  can 
be  identified.  A  nondominated  feasible  decision  is  one  which  does  not  permit  any 
avoidable  idleness  of  the  resources  available  within  a  particular  manpower  structure. 
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Table  2.3:  Manpower  Structure  Alternatives 


9 

h, 

Xi 

X2 

2:3 

Xa 

*5 

Policies 

Iterations 

E{nu) 

R 

1 

1 

2 

1 

2 

0 

0 

4 

1 

.8080 

4.848 

2 

1 

1 

2 

2 

0 

0 

2 

1 

.8159 

4.895 

3 

2 

2 

0 

0 

2 

0 

24 

2 

.7900 

4.740 

4 

2 

1 

0 

0 

3 

0 

12 

2 

.8103 

4.862 

5 

3 

0 

0 

0 

0 

3 

36 

2 

.8409 

5.045 

For  example,  consider  structure  5,  which  provides  three  mechanics  qualified  in  aH 
types  of  maintenance.  Suppose  the  system  is  in  state  15,  so  that  two  aircraft  require 
airframe  and  engine  maintenance.  The  set  D5 15  contains  two  nondominated  feasible 
decisions.  One  option  would  be  to  assign  all  the  mechanics  to  one  aircraft,  with  two 
mechanics  performing  the  engine  task  and  one  mechanic  performing  the  aiirframe 
task  (si5i2  =  2,S]5  14  =  1).  Alternatively,  two  mechanics  could  be  assigned  to 
perform  the  airframe  tasks  on  both  aircraft,  with  one  mechanic  unavoidably  idle 
(^15 12  =  0,Si5n  =  2).  As  indicated  in  Table  2.4,  these  options  correspond  to 
decisions  5  and  7,  so  D515  =  {5,7}.  The  nondominated  feasible  decisions  for  every 
manpower  structure  and  system  state  are  displayed  in  Table  2.5. 

Once  all  nondominated  feasible  decisions  are  identified,  an  efficient  Markov 
decision  algorithm  can  be  employed.  It  can  first  be  noted  that  the  underlying 
Markov  process  is  completely  ergodic,  meaning  that  the  steady-state  probabilities 
for  the  system  states  are  independent  of  starting  conditions  since  all  states  are 
accessible  from  each  other.  Under  this  condition,  the  following  “policy  iteration” 
algorithm  can  be  used  to  converge  on  an  optimal  assignment  policy  for  each  man- 


Table  2.4:  Possible  Manpower  Assignment  Decisions 


k 

Manpower  Assignment 

1 

no  maintenance  reqd. 

1 

521  =  1 

1 

•*32  =  1 

1 

*42  =  1 

1 

*53  =  2,  554  =  1 

2 

*33  =  2,  554  =  0 

3 

*53  =  0,  554  =  1 

k  Manpower  Assignment 


B 


6  I  1  562 

2  562 


1  •S73 

2  573 

3  573 


1  Ss 

8  2  58 

3 


1  «95 

2  595 

3  595 

4  595 

5  595 

6  *95 

7  595 


10  1  5 

2  5 


1  5 

11  2  5 

3  5 


=  2 
=  1 


=  1,  576  =  1 
=  1,  576  =  0 
=  0  ,  576  =  1 


=  ll  586  =  2 
=  1?  «86  =  0 
=  0?  ^86  =  2 


=  1,  597  =  2,  598  =  1 

=  1,  597  =  2,  598  =  0 
=  1)  597  =  0,  598  =  1 
=  1,  597  =  0,  598  =  0 
=  0,  597  =  2  ,  598  =  1 
=  0,  597  =  2,  598  =  0 
=  0,  597  =  0,  598  =  1 


7  =  2 
7  =  1 


2,  5n  8  =  1 

1?  5]1  8  =  0 
0-  ^118  =  1 


1  >*129 

2  5i29 

3  5]29 
12  4  5]29 

5  5i29 

6  5i29 

7  5i29 


3  1  5i38 

2  .8,38 


1  -*1-19 

2  5h9 

3  5i4  9 
4  4  5j49 


^149  ■ 


1  5i512 

2  5i512 

3  5i512 
15  4  5]512 

5  5i512 

6  515  12 

7  «I5]2 

8  Sl512 


=  1?  ^12 10  = 
=  li  >*12  10  = 
=  1,  5]2 10  = 

=  I5  >*12  10  = 

=  0,  5i2 10  = 
=  0,  5i2 10  = 
=  0,  5i2  10  = 


=  4 
=  2 


=  2,  5 
=  2,  5 
=  2,  5 
=  2,  5 
=  0,  5 
=  0,  5 
=  0,  5 


=  4,  5i5  14 
=  4,  5i5i4 

=  4,  5]5h 

=  2,  5i5  14 
=  2,  5i5  14 
=  2,  5i5  14 
=  0i  ^1514 
=  0>  ^15  14 


2,  5i2  11  =  1 
2,  5i2  11=0 
0,  5i2  11  =  1 
0,  5i2  11=0 
2,  5i2  11  =  1 
2,  5i2  11=0 
0?  >*1211  =  1 
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power  structure  g  (see  Appendix  A): 

1.  Initialization.  For  each  state  i  and  each  decision  k  E  Dg,,  compute  transition 
rate  values  as  follows: 

7(u)  =  0 

< 

Select  an  initial  assignment  policy  (a  decision  k'  for  each  state  i). 

2.  Value  Determination.  Let  v,  be  the  relative  value  of  occupying  a  particular 
state  i  under  the  current  policy.  Use  r*’  for  the  current  policy  to  solve  the  set 
of  equations 

f;(no)  =  nj,'^  +  t  =  1, 2, . . . ,  7  (2.20) 

j=i 

for  all  relative  values  v,  and  unknown  £^(nu)  by  arbitrarily  setting  v/  to  zero. 

3.  Policy  Improvement.  For  each  state  i,  find  the  decision  k"  E  Dg,  that  maxi¬ 
mizes  the  expression  r^jVj  using  the  relative  values  of  the  current  pobc}'. 
If  k"  is  unchanged  from  k'  for  all  states,  stop  with  the  optimal  policy.  Other¬ 
wise,  decision  k"  becomes  the  new  current  decision,  so  each  r*'  is  set  equal  to 
T^j  .  Return  to  Step  2. 

The  number  of  iterations  required  by  the  algorithm  can  be  minimized  by  choosing  a 
“good”  initial  policy.  A  suitable  initial  policy  is  fortunately  represented  by  a  set  of 
“greedy”  assignment  decisions;  that  is,  decisions  which  always  favor  the  movement 
of  machines  which  are  closest  to  an  operating  condition.  Such  a  policy  results  Rom 
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a  left-to- right  assignment  of  manpower  to  tasks  as  arrayed  in  Table  2.1.  Transitions 
could  be  similarly  arrayed  for  a  system  of  any  size.  Table  2.3  presents  the  number  of 
iterations  required  to  converge  on  optimal  assignment  policies  when  this  initializa¬ 
tion  rule  is  applied  for  the  “Mike’s  Flying  Club”  example.  Like  the  enumeration  of 
a  particular  policy,  each  iteration  of  the  policy  iteration  algorithm  involves  the  solu¬ 
tion  of  I  simultaneous  equations.  Thus,  the  algorithm  requires  about  8/78  =  10.3% 
of  the  computational  effort  required  by  exhaustive  enumeration. 

Application  of  the  policy  iteration  algorithm  to  the  example  problem  yields 
the  results  shown  in  the  last  two  columns  of  Table  2.3.  The  table  displays  the 
expected  number  of  operating  aircraft  and  sortie  generation  rate  achievable  with 
each  manpower  structure.  It  is  apparent  that  Mike  could  improve  the  performance 
of  his  enterprise  by  replacing  his  five  specialized  mechanics  with  three  more  costly 
mechanics  qualified  in  all  types  of  maintenance  {g  =  5).  This  alternative  would 
produce  a  sortie  generation  rate  of  5.045,  which  is  higher  than  the  maximum  rate 
of  4.848  achievable  with  the  current  manpower  structure  {g  =  1).  The  ontimal 
manpower  structure  represents  a  fully  “generalized”  speci2Jization  strategy. 

2.3  Linear  Programming  Model 

While  the  policy  iteration  algorithm  eliminates  the  need  to  enumerate  solu¬ 
tions  for  every  assignment  policy,  it  still  requires  a  complete  solution  for  each  man¬ 
power  structure.  Assignment  decisions  and  manpower  structures  can  be  considered 
simultaneously  through  a  linear  programming  formulation  of  the  decision  model. 
Linear  programming  (LP)  provides  an  established  method  for  solving  Markov  de¬ 
cision  problems  with  finite  state  spaces  [27],  and  the  general  approach  can  be  effi¬ 
ciently  applied  to  the  manpower  specialization  problem. 
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2.3.1  Concurrent  Approach 

All  candidate  manpower  structures  and  assignment  policies  can  be  eval¬ 
uated  concurrently  in  a  single  linear  program.  First,  decision  variables  Pgik  are 
defined  as  the  joint  probabilities  of  employing  structure  g  G  {1,2, . . . ,  (7},  finding 
the  system  in  state  i  6  {1, 2, . . . , /},  and  selecting  assignment  decision  k  G  Dg,.  An 
LP  model  can  then  be  written  as  follows: 

Maximize 

E(na)  =  E  '“i'Wt  (2-21) 

g=l  ,=1  k€Dg, 

subject  to 


P9‘>  ' 

[  J  j=\keDgj  1 

fl’K(jl)  J 

^  =  1,2,  ...,G  i  =  1, 2, — 1 

(2.22) 

53  53  53  Pgik  =  1 

g=l  ,=1  k^Dg, 

(2.23) 

Pgik>0  g=l,2,...,G  t=l,2,...,/  k  ^  Dg, 

(2.24) 

The  constraints  represented  by  Equations  2.22  ensure  that  global  balance  for  the 
queuing  system  is  satisfied.  Equation  2.23  is  a  normalization  constraint  for  conser¬ 
vation  of  probability.  Excluding  the  non-negativity  restrictions  imposed  by  Equa¬ 
tions  2.24,  the  formulation  will  generate  G(I  —  1)  -I-  1  constraints.  The  number  of 
decision  variables  pgn  will  depend  on  the  total  number  of  nondominated  feasible 
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assignment  decisions  for  all  system  states  and  manpower  structures.  For  exam¬ 
ple,  formulation  of  the  “Mike’s  Flying  Club”  problem  requires  71  constraints  and 
95  decision  variables.  When  this  model  is  implemented  on  a  microcomputer  using 
a  standard  software  package  [6],  a  solution  is  obtained  in  only  a  few  seconds  of 
processing  time.  The  following  decision  variables  have  non-zero  optimal  values: 


Psn  —  .1803 
Ps2i  =  .1803 
P531  =  .1564 
P541  =  -0951 
P553  =  -04  84 


P562  =  .0445 
P573  =  .06  78 
P583  =  .0491 
P595  =  .03  42 
Po  10  2  =  .0219 


P5  1]3  =  -0399 
Pd  12  5  =  .0342 
Ps  13 1  =  .0263 
P5  14  5  =  -0168 
P5  15  3  =  -0047 


The  example  results  provide  a  useful  vehicle  for  examining  some  important 
characteristics  of  the  LP  solution.  First,  it  should  be  noted  that  there  is  exactly  one 
non-zero  decision  variable  for  each  system  state.  The  non-zero  variables  therefore 
represent  the  steady-state  probabilities  tt,.  They  also  identify  the  optimal  assign¬ 
ment  decisions.  While  the  model  permits  randomization  between  decisions,  the 
optimal  solution  will  force  all  probability  for  a  given  state  to  the  decision  which 
produces  the  largest  rate  of  increase  in  the  objective  function.  Similarly,  only  one 
manpower  structure  is  represented  in  the  final  solution  {g  =  5),  even  though  the 
model  permits  randomization  between  structures.  The  global  balance  equations 
ensure  that  all  state  probabilities  for  a  given  structure  and  policy  maintain  con¬ 
stant  relative  proportionate  values.  Therefore,  the  optimal  solution  will  force  all 
probability  into  the  best  structure. 

The  concurrent  LP  provides  a  compact  representation  of  any  manpower 
specialization  problem.  The  model  and  supporting  data  can  be  easily  formatted  in 
a  generalized  algebraic  language  for  computer  implementation,  so  the  formulation 
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can  be  useful  for  small  problems.  When  solved  using  a  standard  simplex  or  revised 
simplex  algorithm,  the  concurrent  formulation  will  normally  require  a  smaller  num¬ 
ber  of  LP  iterations  than  would  be  required  if  the  problem  were  decomposed  to  find 
separate  optimal  solutions  for  each  manpower  structure.  However,  each  LP  itera¬ 
tion  for  the  concurrent  formulation  will  require  a  much  larger  number  of  arithmetic 
operations.  Thus,  the  concurrent  formulation  has  limited  practical  value  since  a 
decomposed  version  of  the  problem  will  generally  require  less  computational  effort. 

Issues  concerning  computational  storage  and  analytical  insight  also  impact 
the  practical  utility  of  the  concurrent  approach.  Every  candidate  manpower  struc¬ 
ture  g  generates  7  —  1  constraints  and  at  least  7—1  decision  variables,  so  a  problem 
involving  a  large  number  of  structures  and  system  states  could  easily  exceed  com¬ 
putational  storage  limits.  Furthermore,  while  solution  of  the  concurrent  LP  yields 
an  optimal  structure  and  associated  system  performance,  it  does  not  provide  di¬ 
rect  insight  on  how  the  optimal  specialization  strategy  compares  with  others  under 
consideration. 

2.3.2  Sequential  Approach 

The  disadvantages  of  the  concurrent  formulation  can  be  overcome  through 
a  new  algorithm  which  sequentially  evaluates  all  candidate  manpower  structures 
but  benefits  from  commonality  of  arithmetic  operations.  First,  a  single  LP  is 
formulated  which  permits  any  assignment  decision  that  is  feasible  for  at  least  one 
structure  under  consideration.  Solution  of  this  LP  wiU  produce  an  upper  bound  on 
the  optimal  system  performance.  Then,  for  each  individual  manpower  structure,  a 
constraint  is  added  which  forces  all  infeasible  decision  variables  to  zero.  Integration 
of  the  new  constraint  into  the  linear  programming  tableau  renders  the  original 
solution  infeasible.  However,  feasibility  can  be  restored  by  performing  dual  simplex 


33 


iterations  [26].  These  iterations  are  terminated  when  either  a  feasible  (optimal) 
solution  is  reached  or  the  objective  value  falls  below  a  known  lower  bound  for  the 
specialization  strategy  employed  by  the  current  structure.  This  lower  bound  is 
established  by  the  best  objective  value  from  previously  evaluated  structures  which 
employ  the  same  specialization  strategy. 

The  algorithm  can  be  accelerated  by  evaluating  all  manpower  structures  in 
order  of  their  estimated  relative  performance  (best  to  worst).  A  good  estimate  for 
relative  performance  of  a  structure  g  can  be  obtained  from  the  parameter 


.k' 


i=l 


7(.j) 


-1 


(2.25) 


where  k'  —  Tmn{k  €  Dgi),  A  structure  with  a  low  value  of  6g  will  tend  to  produce 
favorable  system  performance  since  maintenance-intensive  states  will  have  relatively 
short  transition  times  (inverses  of  transition  rates).  The  ranking  derived  from  6g  is 
only  an  estimate  since  “greedy”  (not  necessarily  optimal)  assignment  decisions  are 
assumed  and  transition  times  are  not  weighted  by  steady-state  probabilities. 

The  sequential  linear  programming  algorithm  can  be  concisely  stated  as 

follows: 


1.  For  each  system  state  i,  define  a  set  of  decisions  Z?,  =  Ug-^Dgi.  Let  p,/t  be  the 
joint  probability  of  finding  the  system  in  state  i  and  selecting  decision  k  E  D,. 
Formulate  and  solve  the  bounding  LP, 


Majdmize 

EM  =  E  E 

i=l 


(2.26) 
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subject  to 

P.k  jnJVo  +  5Z[-^]/i7(.i)  1  =  S  II  inu\.^(j.)Po  +  [-^]/i^0,)| 

Arto,  i  J=1  j  J=lJt€Z)j  I  J 

i=  1,2,...,/- 1  (2.27) 

E  E  P.^  =  1  (2-28) 

i=l  *:€D, 

P.A>0  i  =  l,2,...,/  ike  A  (2.29) 

Store  the  linear  programming  tableau. 

2.  For  each  specialization  strategy  h  €  {1,2,...,^},  define  incumbent  optimal 
performance  values  £(no)/i  with  initial  values  of  zero.  Reorder  all  manpower 
structures  g  E  {1,2,  ...,G}  such  that  <  ^2  <  •  •  •  <  where  each  6g  is 
defined  by  Equation  2.25.  For  g  =  1, . . . ,  G, 

(a)  Define  a  set  of  infeasible  decisions  Dgi  =  {k  :  k  £  Di,k  0  Dgi}.  Modify 
the  LP  tableau  from  Step  1  by  adding  the  constraint 

E  E  =  0  (2-30) 

■=1 

(b)  If  the  addition  of  Equation  2.30  results  in  an  infeasible  current  solution, 
perform  dual  simplex  iterations  [2,  p.  182]  until  E(no)  <  E(no)/,j  or  E(no) 
is  optimal.  If  E{no)  is  optimal,  then  let  E{no)hg  =  E{no). 
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3.  Identify  an  optimal  manpower  manpower  specialization  strategy  h'  such  that 
E{nij)h'  >  E{n(i)h  V  h  G  {1, 2, . . . ,  H}. 

The  sequential  algorithm  minimizes  computational  effort  by  taking  advan¬ 
tage  of  state  relationships  which  are  common  to  different  manpower  structures  and 
assignment  decisions.  For  the  “Mike’s  Flying  Club”  example,  Step  1  of  the  al¬ 
gorithm  requires  15  LP  iterations  to  establish  the  initial  bounding  solution,  and 
Step  2  requires  3-f-l-|-3-f3-|-l  =  11  iterations  to  evaluate  the  individual  man¬ 
power  structures.  The  26  total  iterations  generate  performance  results  for  all  three 
specialization  strategies  under  consideration. 

Two  variations  of  the  sequential  algorithm  may  be  useful  for  particular  ap¬ 
plications.  In  some  cases,  it  may  be  desirable  to  gain  additional  computational 
speed  in  exchange  for  less  analytical  insight.  This  is  accomplished  by  modifying 
Step  2  so  that  all  incumbent  optimal  objective  values  E{n{y)h  are  replaced  by  a  sin¬ 
gle  global  value  £?(no)opf  For  each  manpower  structure  p,  dual  simplex  iterations 
are  performed  until  £?(no)  <  i^(Tio)opt  or  E{no)  is  optimal.  If  E{no)  is  optimal,  it 
becomes  the  new  £^(no)opf  The  modified  algorithm  identifies  only  a  global  opti¬ 
mal  manpower  structure  rather  than  an  optimal  structure  for  each  specialization 
strategy.  However,  execution  of  this  “speed”  version  of  the  algorithm  can  require 
a  significantly  smaller  number  of  LP  iterations. 

Another  variation  of  the  algorithm  can  provide  full  analytical  insight  by 
completely  evaluating  all  manpower  structures.  For  every  structure,  dual  simplex 
iterations  are  continued  until  an  optimal  objective  value  is  achieved.  No  benefit  is 
derived  by  evaluating  structures  in  a  particular  order,  so  the  reordering  procedure 
can  be  deleted.  However,  more  LP  iterations  are  required  by  this  “insight”  version 
of  the  algorithm. 
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All  three  versions  of  the  sequential  linear  programming  algorithm  (normal, 
speed,  and  insight)  are  implemented  in  the  computer  program  described  in  Ap¬ 
pendix  B.  Table  2.6  displays  computational  performance  results  (number  of  LP 
iterations)  for  each  version  when  applied  to  three  example  problems  of  different 
sizes.  Each  iteration  count  includes  all  pivots  to  establish  an  initial  “greedy”  solu¬ 
tion  to  the  bounding  LP,  all  primal  pivots  to  optimize  the  bounding  solution,  and 
aU  dual  pivots  required  by  Step  2  of  the  algorithm.  The  iteration  counts  are  com¬ 
pared  against  a  baseline  number  of  iterations  which  are  required  if  separate  linear 
programs  are  solved  for  each  manpower  structure.  It  is  noteworthy  that  the  normal 
version  of  the  algorithm  requires  28.1-33.8%  of  the  baseline  computational  effort, 
and  the  speed  version  requires  13.4-27.9%  of  the  baseline  effort.  The  LP  iteration 
counts  correlate  approximately  with  comparative  run  times  for  each  problem  size. 


Table  2.6:  Computational  Performance  of  Sequential  Algorithms 


Problem 

B 

G 

Number  of  LP  Iterations 

Insight 

Normal 

Speed 

Mike’s  Flying  Club  {N  =  2,  C  =  100) 

15 

5 

79 

31 

26 

22 

Mike’s  Flying  Club  {N  =  4,  C  =  200) 

mm 

8 

613 

291 

172 

122 

YF-XX  (see  Chapter  3) 

455 

11 

5618 

1110^9 

1897 

774 

2.4  Network  Reduction  Procedure 

Even  when  manpower  structures  are  considered  sequentially,  problem  di¬ 
mensionality  can  still  be  an  important  issue.  For  example,  consider  a  system  with 
three  machines  and  six  task  types,  where  all  tasks  can  be  accomplished  simul¬ 
taneously.  The  corresponding  queuing  network  will  have  Z  —  =  255 
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maintenance  conditions,  so  the  system  will  have  I  =  {^^55°)  —  2.829  x  10®  states. 
Fortunately,  the  routing  probabilities  and  average  service  times  associated  with 
many  of  the  maintenance  conditions  could  be  very  small,  so  a  very  good  approxi¬ 
mate  solution  could  be  obtained  by  “collapsing”  these  conditions  into  others.  The 
notion  of  flow  equivalence  can  be  applied  to  form  a  new  reduced  network  that  wiU 
approximate  the  behavior  of  the  original  network. 

Like  all  previous  concepts,  the  network  reduction  approach  can  be  demon¬ 
strated  using  “Mike’s  Flying  Club.”  Suppose  computational  facilities  are  so  limited 
that  any  system  with  more  than  ten  states  is  intractable.  It  is  therefore  necessary 
to  develop  a  new  network  that  will  approximate  the  original  example  network  while 
reducing  the  total  number  of  system  states.  This  can  be  done  by  eliminatng  sta¬ 
tion  4  and  adjusting  the  characteristics  of  the  stations  which  can  be  entered  directly 
from  station  4  (stations  2  and  3).  The  adjusted  characteristics  include  the  rout¬ 
ing  probabilities  and  service  rates.  The  resulting  reduced  network  will  have  three 
maintenance  stations  and  (3)  =  10  states. 

Routing  probability  can  be  conserved  by  adjusting  92  ajid  to  new  values 
92  and  q'^  as  follows: 


92 


=  92  +  94  (—^) 

V^2  +  /*3/ 

=  .1404  + (.1454)  (:^) 


=  .2373 


93  —  93  + 


\/i2  + 


=  .1880 +  (.1454)  (33!^)  = 


.2364 


This  operation  distributes  94  to  succeeding  stations  in  proportion  to  the  relative 
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probabilities  of  machine  entry  from  station  4.  It  is  now  necessary  to  decrease  the 
service  rates  at  stations  2  and  3  to  reflect  the  longer  service  times  for  machines 
requiring  both  types  of  maintenance; 


/ 

/ij 


f  (i)  +  (i  ^ 


-1 


+  9-1  (^) 


-  I® 


^3 


+  M3, 


'  1  1 

—  +  - 

,M2  M2  +  M3, 


.2373 

1.0785 


=  .2200 


-1 


M3 


^  (i)  (;^)  (ot  + 


-1 


93  1 93 

.2364 

.5376 


:a-( 


M2 


1 


+ 


M2-\-M3/  \M3  M2  +  M3, 


=  .4400 


The  reduction  in  service  rates  compensates  for  the  deletion  of  station  4. 

This  example  reduction  can  be  expanded  into  a  general  procedure.  First, 
let  /'  <  /  be  a  maximum  number  of  system  states  imposed  by  computational 
limitations.  The  maximum  number  of  reduced  network  stations  can  therefore  be 
represented  as  Z'  <  Z,  where  )  <  /'  <  If  all  maintenance  stations 

are  ordered  by  increasing  numbers  of  tasks  in  their  pending  task  sets,  a  network 
reduction  procedure  can  be  constructed  as  follows: 

1.  Deflne  station-dependent  task  rates  Mm^  f<>^  each  station  and  applicable  task. 
Initially,  let  Mm^  =  Mm  for  each  z  €  {1, 2, . . . ,  Z}  and  each  eligible  task  m  ^  E^. 
Define  z~  as  the  current  station  to  be  removed  from  the  network.  Initially,  let 
2-  =  Z. 
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2.  Let  5j-  be  the  set  of  stations  which  can  be  directly  entered  from  station  z~ . 
For  each  station  z'^  G  S^- , 

(a)  Identify  the  task  e  G  E,-  which  causes  machine  movement  to  station  z"*" 
on  completion. 

(b)  Compute  a  new  routing  probability  as 


9z+  — 


l-imeE  _ 


(2.31) 


(c)  For  each  task  m  G  ,  compute  a  new  service  rate  eis 


/  (j"*' ) 


“fn 


+  qz- 


Mi- 


/ 

/ 

V 

> 

-  M»)  ^ 

M») 

1 


(2.32) 


1" 


3.  For  each  z"*"  in  5j-,  let  9^+  =  and  V  m  G  E^-t^.  Let 

z~  =  z~  —  1.  If  z“  =  Z',  then  stop  with  the  final  reduced  network.  Otherwise, 
return  to  Step  2. 


This  procedure  will  transfer  maintenance  time  at  stations  with  many  pend¬ 
ing  tasks  to  stations  with  fewer  pending  tasks.  Consequently,  the  reduced  network 
will  not  capture  the  queiung  delays  caused  by  resource  conflicts  between  the  eligi¬ 
ble  tasks  at  each  eliminated  station.  The  approximation  will  therefore  overestimate 
system  performance.  However,  an  approximate  flow  equivalent  network  can  pro¬ 
vide  very  accurate  results  unless  multiple  task  eligibilities  are  very  common  and 
resource  availability  is  very  low.  When  solving  the  example  problem  using  the  re¬ 
duced  network,  an  optimal  sortie  generation  rate  of  5.130  sorties  per  aircraft  per 
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day  is  obtained.  This  represents  an  error  of  only  1.75%,  even  though  a  sizable 
portion  ^1454)  of  the  network  routing  probability  is  redistributed.  More  impor¬ 
tantly,  the  optimal  manpower  structure  is  unchanged  from  the  full  network  solution, 
x  =  (0,0, 0,0, 3). 


Chapter  3 


Application 

In  the  previous  chapter,  a  simple  example  was  used  to  develop  a  queuing 
network  approach  to  maintenance  manpower  specialization.  In  this  chapter,  the 
method  is  applied  to  a  larger  problem  which  is  more  representative  of  potential 
“real  world”  applications.  The  problem  objective  is  to  derive  an  optimal  manpower 
structure  for  dispersed  operation  of  a  notional  military  aircraft. 

3.1  The  YF-XX  Aircraft  Maintenance  Problem 

The  YF-XX  is  a  notional  tactical  fighter  in  the  prototype  phase  of  weapon 
system  development.  Since  the  earliest  stages  of  the  system  design,  significant  effort 
has  been  focused  on  reliability  and  maintainability  characteristics.  Consequently, 
substantial  analysis  has  been  conducted  to  estimate  subsystem  failure  rates,  repair 
times,  and  task  personnel  requirements.  Since  extensive  test  data  has  not  yet 
been  accumulated,  precise  forms  for  probability  distributions  of  failure  and  repair 
times  are  not  known.  However,  efforts  to  eliminate  known  failure  mechanisms 
suggest  that  subsystem  malfunctions  can  be  regarded  as  “random”  and  can  thus  be 
accurately  represented  by  exponentially  distributed  f^ure  times.  The  exponentiad 
distribution  is  also  postulated  as  a  suitable  model  for  maintenance  task  time  at  the 
major  subsystem  level.  The  aggregate  maintenance  data  for  the  YF-XX  are  shown 
in  Table  3.1. 

The  operational  concept  for  the  YF-XX  involves  a  wide  variety  of  mission 
profiles  with  an  average  duration  of  1.6  hours  (/xq  =  .625).  Tum-around  mainte- 
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Table  3.1;  Task  Data  for  YF-XX  Tactical  Fighter 


m 

Task  Description 

1 

Munitions  Upload 

- 

2.400 

3 

2 

Aircraft  Turn-Around 

- 

2.000 

1 

3 

Avionics  Repair 

.0713 

.4191 

1 

4 

General  Aircraft  Repair 

.0269 

.2390 

2 

5 

Engine  Repair 

.0250 

.2301 

2 

6 

Electrical  Subsystem  Repair 

.0232 

.2107 

1 

7 

Pneudraulic  Subsystem  Repair 

.0144 

.4223 

2 

8 

Fuel  Subsystem  Repair 

.0103 

.1812 

2 

9 

Armament  Subsystem  Repair 

.0091 

.3200 

2 

nance  and  munitions  upload  are  accomplished  between  all  sorties.  AU  unsched¬ 
uled  maintenance  generated  by  subsystem  malfunctions  must  be  performed  prior 
to  aircraft  turn-around,  and  turn-around  must  be  completed  before  munitions  up¬ 
load.  Additionally,  all  electrical  subsystem  repair  must  precede  any  work  on  the 
aircraft  avionics  (communications,  navigation,  electronic  counter-measures,  etc.). 
Figure  3.1  displays  a  reduced  network  which  captures  these  requirements.  The  fig¬ 
ure  also  displays  the  adjusted  routing  probabilities  associated  with  each  network 
station.  While  the  unreduced  network  included  129  maintenance  stations,  less 
than  2.5%  of  the  total  routing  probability  had  to  be  redistributed  among  stations  3 
through  12  using  the  procedure  described  in  Section  2.4. 

3.2  Maintenance  Cost  Model 

The  basic  maintenance  concept  for  the  YF-XX  assigns  the  nine  task  types 
to  seven  manpower  specialties.  However,  wartime  operation  plans  dictate  that  a 
significant  number  of  aircraft  will  operate  in  groups  of  three  at  dispersed  locations. 
It  is  anticipated  that  some  maintenance  specialty  merger  will  be  required  for  dis- 
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persed  forces  in  order  to  satisfy  manpower  cost  constraints.  Table  3.2  lists  thirteen 
types  of  specialties  which  appear  in  various  consolidation  schemes  under  consid¬ 
eration.  The  table  displays  the  task  qualifications  for  each  specialty  and  provides 
additional  information  that  can  be  used  to  develop  cost  estimates.  The  estimated 
length  (in  days)  of  the  initial  training  period  for  each  type  of  specialty  is  indicated 
as  Ly.  All  training  generates  a  direct  cost  of  $200  per  person  per  day.  Also  shown 
is  the  expected  monetary  bonus  Vy  (in  dollars)  which  would  be  paid  to  each  type 
of  specialist  if  he  were  to  reenlist  after  a  four-year  term  of  service.  This  value  is 
largely  determined  by  the  demand  for  personnel  with  similar  skills  in  the  civilian 
economy. 

An  expected  profile  of  the  entire  maintenance  workforce  is  presented  in 
Figure  3.2.  This  profile  is  adapted  from  a  demonstrative  model  used  in  the  Air  Force 
SUMMA  project  [28,  p.  10].  The  model  asserts  that  about  60%  of  the  workforce  will 


44 


Table  3.2:  Maintenance  Specialties  for  YF-XX 


y 

Qy 

Ly 

^y 

1 

{1,9} 

90 

- 

20,934 

2 

{2,4} 

100 

5,000 

21,637 

3 

{3} 

120 

10,000 

22,679 

4 

{5} 

120 

5,000 

22,302 

5 

{6} 

95 

10,000 

21,851 

6 

{7} 

90 

- 

20,934 

7 

{8} 

90 

- 

20,934 

8 

{3,6} 

0 

00 

10,000 

24,799 

9 

{5,8} 

150 

5,000 

23,338 

10 

{5,7,8} 

200 

5,000 

25,174 

11 

{1,2, 4,9} 

150 

5,000 

23,338 

12 

{1,2,4,5,7,8,9} 

270 

12,000 

28,531 

13 

{1,2,3,4,5,6,7,8,9} 

360 

15,000 

32,924 

depart  after  a  single  four-year  term  of  service,  and  40%  of  the  remaining  workforce 
will  depart  after  a  second  term  of  service.  Bonuses  are  to  be  paid  for  second 
and  third  term  reenlistments  in  order  to  restrict  workforce  attrition  to  these  levels 
for  all  specialties.  It  is  further  postulated  that  all  personnel  who  have  completed 
three  terms  of  service  will  remain  until  retirement  after  a  20-year  career.  These 
parameters  establish  portions  of  the  workforce  Wr  for  each  term  t  6  {1,2, ...  ,5}. 
The  average  annual  cost  of  pay  and  benefits  Ur  for  personnel  in  each  term  is  also 
shown.  Pay  and  benefits  do  not  vary  with  specialty. 

Training  costs,  retention  costs,  and  the  information  presented  in  Figure  3.2 
can  be  aggregated  to  produce  total  annual  costs  for  each  specialty  as  follows: 


=  Wi 


-I- 


(  W(365)  \(‘^r  +u\ 

l,(4)(365) -lJ  \  4  + 

W2  ^2  +  4"  ^3  (^3  4"  +  W^Ua  +  W3U3 


(3.1) 
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Workforce 

Density 


Years  of  Service 


Figure  3.2:  Profile  of  YF-XX  Manpower  Force 


The  aggregate  cost  parameters  listed  in  the  last  column  of  Table  3.2  are  computed 
using  this  formula.  Note  the  the  first  term  of  the  equation  includes  an  expression 
which  accounts  for  the  opportunity  cost  of  training  time.  Clearly,  training  and 
retention  factors  cause  a  significant  increase  in  the  unit  cost  of  maintenance  per¬ 
sonnel  as  specialization  decreases.  Table  3.2  indicates  that  the  annual  cost  of  a 
maintenance  technician  with  a  full  range  of  skills  is  roughly  50%  more  than  that  of 
a  typical  technician  with  a  highly  specialized  skill. 

3.3  Manpower  Optimization 

The  annual  maintenance  manpower  budget  for  the  YF-XX  is  constrained  to 
a  level  of  $120,000  per  aircraft.  Since  wartime  plans  assume  that  each  maintenance 
technician  will  be  on  duty  for  a  12-hour  shift  daily,  the  total  expenditure  limit 
per  shift  at  each  dispersal  base  can  be  calculated  as  C7  =  (3)(12/24)(120,000)  = 
180,000.  As  indicated  in  Table  3.3,  this  constraint  translates  to  a  total  of  11  can¬ 
didate  manpower  structures.  Note  that  the  expenditure  limit  will  not  support  a 
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fully  specialized  structure.  Several  specialty  types  do  not  appear  in  any  of  the 
acceptable  alternatives. 


Table  3.3:  Manpower  Structures  for  YF-XX 


*  local  optimum  for  specialization  strategy 
**  global  optimum 


Table  3.3  also  displays  the  performance  results  obtained  when  queuing  net¬ 
work  analysis  is  applied.  The  reduced  network  shown  in  Figure  3.1  generates 
=  455  system  states,  so  the  problem  is  quite  tractable  when  solved  using 
the  sequential  linear  programming  algorithm.  Results  are  shown  for  each  man¬ 
power  structure,  indicating  a  steady  improvement  in  sortie  generation  capability 
as  specialization  decreases.  Dispersal  bases  will  operate  most  effectively  if  each 
maintenance  shift  is  manned  by  five  personnel  that  are  qualified  on  all  aircraft  sys¬ 
tems.  Structure  8  also  produces  good  results  and  might  be  favored  in  the  context 
of  broader  manpower  issues.  This  structure  employs  one  avionics/electrical  system 
specialist,  and  five  technicians  qualified  in  all  other  types  of  maintenance. 


Chapter  4 


Extensions 

The  general  model  presented  in  this  thesis  can  be  adapted  to  capture  a 
variety  of  maintenance  concepts  which  might  arise  in  specific  applications.  In  this 
chapter,  the  versatility  of  the  model  is  demonstrated  by  extending  it  to  address 
some  issues  which  might  apply  to  the  aircraft  maintenance  problem  represented  by 
the  “Mike’s  Flying  Club”  example. 

4.1  Cross-Training 

It  has  been  shown  that  a  queuing  network  model  can  be  used  to  optimize 
the  basic  structure  of  a  maintenance  manpower  force.  This  same  approach  can  be 
extended  to  determine  optimal  “cross-training”  for  maintenance  specialists  so  that 
they  possess  secondary  skills. 

Suppose  the  “Mike’s  Flying  Club”  enterprise  currently  employs  a  fully  spe¬ 
cialized  maintenance  force  consisting  of  two  turn-around  mechanics,  one  airframe 
mechanic,  and  two  engine  mechanics.  Maintenance  specialists  with  broader  skills 
are  not  immediately  available,  but  Mike  has  the  opportunity  to  provide  some  em¬ 
ployees  (current  or  future)  with  additional  training.  Specifically,  he  can  qualify  one 
or  more  airframe  mechanics  to  assist  an  engine  specialist  with  engine  repairs.  The 
cost  of  this  training  would  translate  into  an  incremental  expense  of  $10  per  hour. 
In  addition,  cross-trained  airframe  mechanics  could  be  further  upgraded  to  perform 
turn-around  functions  at  an  incremental  cost  of  only  S3  per  hour. 

The  cross-training  options  selected  for  this  example  result  in  a  set  of  po- 
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tential  maintenance  specialties  with  unit  costs  and  task  qualifications  which  are 
identical  to  those  used  in  earlier  analysis  (see  Table  2.2).  However,  the  concept 
of  cross-training  produces  a  very  different  set  of  acceptable  manpower  structures. 
The  earlier  requirement  that  a  maintenance  task  can  be  allocated  to  only  one  type 
of  specialist  is  now  eliminated.  Instead,  an  acceptable  manpower  structure  must 
satisfy  the  following  constraints: 

1.  Total  hourly  manpower  expenditure  must  not  exceed  the  fixed  budget;  i.e., 
Sy=l 

2.  Sufficient  mechanics  must  be  qualified  to  perform  or  assist  in  each  type  of  task; 
i.e.,  ^y,m^Qy  Xy  >  a^n  Tn  E  {1,2,3}. 

3.  Sufficient  primary-skilled  mechanics  must  be  available  for  each  type  of  task 
requiring  multiple  personnel;  i.e.,  «3  >  1  (at  least  one  engine  mechanic  must 
be  available). 

4.  AU  skills  must  have  a  potential  for  utilization;  i.e., 

12y:meQy  ®V  —  ^rnN  Vm  G  {1,2,3}. 

These  constraints,  in  general  form,  will  apply  to  a  problem  of  arbitrary  size.  The 
fact  that  cross-trained  specialists  may  not  be  as  effective  in  their  secondary  skills  as 
they  are  in  their  primary  skills  can  be  reflected  in  lower  task  completion  rates.  For 
this  example,  assume  that  an  engine  mechanic  and  cross-trained  airframe  mechanic 
can  accomplish  engine  repairs  at  a  slightly  degraded  rate  of  .45  tasks  per  hour.  Sim¬ 
ilarly,  assume  that  a  cross-trained  airframe  mechanic  can  accomplish  turn-around 
tasks  at  an  average  rate  of  .9  per  hour.  Table  4.1  summarizes  the  rate  capabilities 
for  each  type  of  specialist  and  lists  all  of  the  candidate  manpower  structures. 
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Table  4.1:  Manpower  Structures  with  Cross-Training 


9 

®2 

13 

®5 

-E(no) 

R 

1 

2 

1 

2 

0 

0 

.8080 

4.848 

2 

1 

2 

2 

0 

0 

.8159 

4.895 

3 

2 

1 

1 

1 

0 

.8322 

4.993 

4 

2 

0 

2 

1 

0 

.8080 

4.848 

5 

1 

1 

1 

0 

1 

.8302 

4.981 

6 

1 

0 

2 

0 

1 

.8062 

4.837 

7 

1 

0 

1 

2 

0 

.7970 

4.782 

8 

1 

0 

1 

1 

1 

.8302 

4.981 

9 

0 

0 

1 

0 

2 

.8065 

4.839 

/^i 

1.0 

0.90 
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0.25 

0.25 
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The  policy  iteration  or  linear  programming  optimization  methods  can  be 
applied  to  solve  a  cross-training  problem  with  minor  adaptation.  The  set  of  man¬ 
power  assignment  decisions  for  each  state  may  expand,  since  some  decisions  may 
need  to  reflect  a  choice  of  which  specialist  type(s)  are  assigned  to  each  transition 
task.  Transition  rates  are  determined  not  only  by  manpower  assignment  variables 
sfj  ,  but  also  by  corresponding  task  completion  rates 

The  last  two  columns  of  Table  4.1  display  the  performance  results  for  the 
“Mike’s  Flying  Club”  problem  with  cross-training.  The  optimal  sortie  generation 
rate  of  4.993  is  achieved  with  manpower  structure  3,  where  x  =  (2, 1, 1, 1,0).  Mike 
should  therefore  dismiss  an  engine  mechanic,  hire  another  airframe  mechanic,  and 
provide  cross-training  for  one  airframe  mechanic  in  engine  repairs.  This  new  man¬ 
power  structure  will  improve  his  average  sortie  generation  rate  by  about  3%  over 
the  original  rate  achieved  with  structure  1. 
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4.2  Dependent  Machine  Operations 

Another  modeling  problem  which  may  be  encountered  in  applications  is  a 
need  to  operate  two  or  more  machines  together  (e.g.,  flights  of  multiple  aircraft). 
Suppose  half  of  the  flights  scheduled  in  the  “Mike’s  Flying  Club”  enterprise  are  ac¬ 
robatic  formation  flights  which  require  two  aircraft.  This  operational  requirement 
can  be  modeled  using  the  network  shown  in  Figure  4.1.  The  network  station  O' 
represents  the  operating  condition  for  single  aircraft  flights,  while  station  0^  repre¬ 
sents  the  operating  condition  for  formation  flights.  Station  0^  is  also  used  to  hold 
a  single  aircraft  that  has  completed  maintenance  and  is  waiting  for  a  companion 
aircraft.  In  general,  a  network  station  0*'  is  required  for  every  group  size  b  E  B, 
where  B  is  the  set  of  all  required  group  sizes  {B  =  {1,2}  for  the  example).  The 
expected  number  of  operating  machines  can  be  computed  as 


0*:  aircraft  holding/operating 
in  formation 


O':  aircraft 
operating 

1:  turn-around 
muntenance 

2:  airframe 
maintenance 

3:  engine 

maintenance 

4;  airframe  and 
maintenance 


Figure  4.1:  Dependent  Aircraft  Operation 
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Table  4.2  describes  the  21  states  for  the  example  system  with  N  =  2  aircraft 
and  illustrates  the  feasible  state  transitions.  Transition  activities  for  the  lower 
right  portion  of  the  table  are  identical  to  those  for  the  example  without  dependent 
aircraft  operations  (see  Table  2.1).  However,  transition  rates  involving  maintenance 
tcLsk  1  now  reflect  the  different  stations  which  an  aircraft  can  enter  from  station  1. 
Define  p*  =  .5  as  the  probability  that  the  next  scheduled  flight  requires  a  single 
aircraft.  Similarly,  define  =  .5  as  the  probability  that  the  next  scheduled  flight 
requires  two  aircraft.  Transition  rates  can  be  expressed  as  for 

(ij)  e  {(8, 7),  (12, 8),  (13, 9),  (14, 10),  (15, 11)}  and  as  for  (i,i)  € 

{(8, 2),  (12, 3),  (13, 4),  (14, 5),  (15, 6)}.  For  the  case  where  an  aircraft  is  holding  for 
a  formation  companion  (state  3),  probabilistic  branching  does  not  apply  since  the 
next  aircraft  leaving  station  1  must  enter  station  0^. 


Table  4.2:  Transition  Activities  (7)  for  Dependent  Operations 
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While  a  single  gaining  station  z'*'(ij)  was  formerly  associated  with  each 
transition,  this  parameter  must  now  be  replaced  by  a  list  of  gaining  stations 
The  routing  probabilities  9^+  must  also  be  displaced  by  new  values  5^+  for  each 
station  0^.  The  computation  of  a  routing  probability  5^+  **'  single  aircraft 

completing  a  sortie  is  unchanged  from  previous  analysis,  but  more  effort  is  required 
to  obtain  probabilities  for  aircraft  completing  a  formation  flight.  Let  T02  be  the 
duration  of  the  flight,  let  TI2  be  the  time  until  an  airframe  malfunction  occurs 
for  the  first  aircraft  in  the  formation,  and  let  TI3  be  the  time  until  an  engine 
malfunction  occurs  for  the  first  aircraft.  Similarly,  let  T22  and  T23  represent  the 
malfunction  times  for  the  second  aircraft.  Routing  probabilities  5^+  can  then  be 
computed  as  follows: 


fJ-o 


=  P{To2  <  min(ri2,ri3,r22,r23)}  =  =  -3571 


-2 

9{1.2} 


=  P{To^  <  min(Tl2,  Tlsj T23),  Tqz  >  T22} 

+  P{To7  <  min(n3, 722,723),  7o2  >  7I2} 

=  ^  +  a7+  2A,  ■ 


9{1.3} 


P{7o2  <  min(7l2,722,723),  Tq/i  >  7I3} 

+  P{Tq7  <  min(7l2, 722, 7I3),  7o2  >  723} 


+  2A2  +  A3 


...))  - 


1553 


9{i.4}  —  ^{^0=  <  niin(7l2, 7I3),  7o2  >  max(722,723)} 


+  P{7o2  <  min(722, 723),  7o2  >  max(7l2,7l3)} 
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=  2( - _  to)  =  .0641 

\/io  +  ^2  +  2A3  2  2  y 

9f2.2}  ~  -^{^0*  2^23),  To?  ^  max(Tl2,  r22)} 


„2 

5{2.3} 


=  P{To?  <  min(Tl2,T23),  To?  >  max{Tl3,  T22)} 


+  P{To?  <  inin(T22,Tl3),  To?  >  max(Tl2,T23)} 


=  2 


Mo 


^Mo  +  -^2  +  -^3 


=.0641 


9(2.4} 


P{To?  <  TI3,  To?  >  max(Tl2,T22,T23)} 

+  P{To?  <  T23,  To?  >  max(Tl2,  TI3,  T22)} 
c^(  Mo  2  2  9(1.3}  9(1.4}  2  9f2.3}\ 

— r~  — 


=  .0500 


9(3.3}  —  <  niin(Tl2,  T22),  To?  >  max{Tl3,  T23)} 


Mo 


Mu  +  2A2 


=  .0431 


9{3.4} 


=  P{To?  <  TI2,  To?  >  max(Tl3,T22,T23)} 


+  P{To?  <  T22,  To?  >  max(Tl2,  TI3,  T23)} 

^  9fl.4}  2  9{2.3}^ 

—  -  9(3,3}  -  —  j 


9  f  Mo  2  9{i.2}  2 


=  .0704 


„2 

9(4.4} 


P{^o2  ^  ni^(Tl2,Tl3,  T22,T23)} 

1  ~  9(1,1}  ~  9(1.2}  ~  9(1,3}  ~  9(1.4} 

■“  9(2.2}  ~  9(2,3}  ~  9(2,4}  ~  9(3.3}  ~  9(3,4}  =  -0532 
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Unlike  previous  versions  of  the  queuing  network  model,  the  dependent  ma¬ 
chine  extension  does  not  exhibit  one-step  behavior.  However,  since  the  policy  itera¬ 
tion  and  linear  programming  approaches  are  based  on  global  balance,  both  solution 
methods  can  still  be  applied.  For  notational  convenience,  rates  are  defined  for 
every  possible  system  transition: 


3  :  z-(ij)  =  0‘' 


=  < 


3  b' eB:  Z^iij)  =  {0^'}, 
<2  -  =  0  V  fc  €  B 


[  [■sfj7®7(.j)]/i7(.j)  otherwise 


An  extended  version  of  the  concurrent  linear  programming  formulation  can  then 
be  written  as  follows: 


Maximize 


c  /  _(') 

g=]  ,=1  k^Dg- b€B 


(4.1) 


subject  to 


Pfl**  a  ’’o  =  12  Psjk 

k€Dg,  J=\  J=1  k^Ogj 


1,2,. ..,G 


i  =  1,2,...,/- 1  (4.2) 


G  / 


1212  12  P9ik  =  1 

9=1  «=>  fceOj, 


(4.3) 


Pffik  S3  S  =  1  —  1,2,...,/  /c£  Dg, 


(4.4) 
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Table  4.3  displays  performance  results  for  all  manpower  structures  in  the 
extended  “Mike’s  Flying  Club”  example.  It  is  assumed  that  average  sortie  dura¬ 
tions  are  identical  for  single  aircraft  flights  and  formation  flights  (/io‘  =  =  .5), 

although  this  is  not  required  in  the  general  model.  Sortie  generation  rates  with  and 
without  dependent  machine  operations  are  shown  comparatively  with  performance 
ranks  in  parenthesis.  Note  that  overall  performance  declines  significantly  with  de¬ 
pendent  machine  operations.  This  loss  can  be  attributed  to  the  idle  time  expended 
by  ready  aircraft  that  must  wait  for  a  formation  companion.  While  the  fully  gen¬ 
eralized  manpower  structure  is  optimal  for  each  case,  the  performance  ranks  for 
the  other  structures  vary.  Dependent  machine  operations  thus  have  a  significant 
impact  on  the  relative  merit  of  the  various  manpower  structures. 


Table  4.3:  Results  with  Dependent  Machine  Operations 


Chapter  5 


Conclusions 

In  a  survey  of  analytical  methods  for  cyclic  queues  and  closed  queuing 
networks,  a  prominent  theoretician  observed  that  “despite  the  extent  of  their  ap¬ 
plications,  the  richness  and  explanatory  power  of  these  methods  have  not  been 
appreciated  or  understood  by  many  practitioners  ”  [23,  p.  605].  This  thesis  has 
demonstrated  that  the  conceptual  framework  of  a  closed  queuing  network  can  offer 
valuable  insight  into  a  “real  world”  problem  which  practitioners  have  tradition¬ 
ally  approached  through  other  methods.  Through  specified  modeling  assumptions 
which  enable  the  use  of  Markov  decision  analysis,  a  queuing  network  model  can  be 
employed  to  determine  an  optimal  level  of  specialization  and  optimal  task  alloca¬ 
tion  for  a  maintenance  manpower  force.  The  method  can  be  specifically  applied 
to  the  problem  of  maximizing  the  operational  effectiveness  of  a  dispersed  unit  of 
military  aircraft  subject  to  a  constraint  on  maintenance  manpower  expenditure. 

Two  important  techniques  can  be  applied  to  solve  a  queuing  network  model 
for  manpower  optimization:  policy  iteration  and  linear  programming.  The  policy 
iteration  approach  allows  an  analyst  to  conjecture  optimal  task  assignments  for 
each  manpower  structure  under  consideration.  If  these  conjectured  assignments 
are  good,  a  policy  iteration  algorithm  will  quickly  converge  to  a  solution.  The 
policy  iteration  approach  is  particularly  useful  for  problems  where  the  number  of 
candidate  manpower  structures  is  relatively  small.  The  main  disadvantage  of  this 
technique  is  that  at  least  one  set  of  simultaneous  equations  must  be  solved  for  every 
manpower  structure  under  consideration. 
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Linear  programming  is  another  powerful  tool  which  can  be  marshalled  to 
find  an  optimal  specialization  strategy.  A  concurrent  LP  formulation  can  theoret¬ 
ically  be  developed  to  simultaneously  consider  all  candidate  manpower  structures 
and  assignment  policies  in  a  single  model.  The  concurrent  model  is  mathemat¬ 
ically  simple,  but  computational  inefficiency  and  high  storage  requirements  limit 
the  utility  of  this  approach  for  larger  problems.  These  shortcomings  can  be  avoided 
through  an  algorithm  which  sequentially  solves  models  for  each  structure  but  elim¬ 
inates  much  of  the  computational  effort  required  for  each  solution.  The  sequen¬ 
tial  approach  also  provides  full  insight  into  the  relative  merits  of  all  specialization 
strategies  under  consideration. 

The  queuing  network  approach  to  manpower  optimization  cannot  com¬ 
pletely  displace  all  other  evaluation  methods.  Some  required  assumptions  may  not 
be  appropriate  for  a  particular  application.  The  extensions  presented  have  demon¬ 
strated  that  creative  modeling  techniques  can  increase  the  flexibility  of  queuing 
network  an«dysis,  but  the  range  of  this  flexibility  is  constrained.  Another  limita¬ 
tion  is  imposed  by  problem  complexity.  The  queuing  network  model  suffers  from 
the  “curse  of  dimensionality”  [3,  p.  323].  As  the  numbers  of  machines  and  task 
types  increase,  the  system  state  space  enlarges  rapidly.  Approximation  methods 
which  compromise  the  “discreteness”  of  machines  in  the  system  are  not  helpful 
because  they  mask  an  important  characteristic  of  the  problem.  The  use  of  approx¬ 
imate  flow  equivalence  to  reduce  network  size  offers  vital  relief,  but  the  practical 
appbcability  of  the  queuing  network  approach  is  still  restricted  to  simple  s}"stems 
or  highly  aggregated  analysis  of  complex  systems. 

Problem  size,  violations  of  assumptions,  or  the  need  to  capture  the  effects  of 
other  dimensions  of  an  overall  logistics  system  may  force  a  primary  reliance  on  sim¬ 
ulation  modeling.  However,  a  queuing  network  model  may  provide  vjiluable  insight 
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when  used  as  an  adjunct  to  simulation.  For  example,  the  analytical  model  can  de¬ 
termine  the  assignment  policies  to  be  employed  by  candidate  manpower  structures 
in  a  simulation  study.  The  simple  “Mike’s  Flying  Club”  problem  would  require 
78  simulation  experiments  to  exhaustively  evaluate  all  policies  and  structures,  and 
multiple  replications  would  be  required  to  place  statistical  confidence  bounds  on 
each  result.  Preliminary  queuing  network  analysis  could  instantly  narrow  the  range 
of  policies  which  merit  investigation.  The  analytical  model  could  also  be  employed 
to  produce  “external  control  variates”  for  increasing  the  computational  efficiency 
of  a  simulation  effort  [25,  p.  359].  Other  adjunct  uses  are  possible,  such  as  using 
the  analytical  model  to  identify  a  starting  point  for  a  simulation  search. 

The  principal  value  of  the  queuing  network  method  fies  is  its  ability  to 
quickly  reveal  the  relationships  between  key  parameters,  evaluate  tradeoffs,  offer 
fundamental  insights,  ajid  answer  basic  questions.  Listed  below  are  some  ques¬ 
tions  which  might  arise  in  the  context  of  manpower  analysis  for  military  aircraft 
maintenance: 

•  How  does  a  change  in  the  manpower  budget  affect  operational  effectiveness? 
Is  the  optimal  specialization  strategy  altered?  How  much  effectiveness  is  lost  if 
the  budget  declines  but  the  specialization  strategy  remains  unchanged?  What 
effects  result  from  changes  in  manpower  training  costs? 

•  What  effect  does  changing  the  reliability  or  maintainability  characteristics  of 
a  particular  aircraft  subsystem  have  on  operational  effectiveness?  What  is  the 
impact  of  such  a  change  on  the  optimal  specialization  strategy?  How  does  the 
cost  of  improving  the  subsystem  compare  with  the  manpower  costs  that  might 
be  saved? 

•  What  is  the  cost  of  dispersed  operations?  How  much  sortie  generation  capa- 
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bility  is  lost  if  aircraft  are  dispersed  without  an  increase  in  total  manpower? 
How  much  of  an  increase  in  manpower  expenditure  is  required  to  maintain  a 
specified  sortie  generation  rate? 

These  questions  are  representative  of  myriad  possible  concerns.  Similar  issues  might 
arise  for  many  other  enterprises  which  employ  maintenance  personnel  of  varying 
skills.  These  issues  can  be  explored  through  the  analytical  method  contributed  by 
this  thesis. 


Appendix  A 


Policy  Iteration  Algorithm 

This  appendix  includes  a  derivation  and  convergence  proof  for  the  policy 
iteration  algorithm  presented  in  Section  2.2.  The  algorithm  can  be  described  as  a 
special  case  of  the  general  policy  iteration  method  developed  by  R.  Howard  [18]. 

Consider  a  machine  maintenance  system  with  Markovian  transition  rates 
r,j  defined  as  follows: 

7(u)  =  0 

0  otherwise 

During  the  time  the  system  is  in  state  i,  it  accumulates  machine  operating  time  at 
a  rate  Let  dt  be  an  infinitesimal  time  interval  such  that  a  system  in  state  i 
will  transition  to  state  j  in  dt  with  probability  r.ydt.  Let  be  the  probability 

that  the  system  will  be  in  state  j  at  time  t  after  the  start  of  the  Markov  process. 
Further,  let  v,{t)  be  the  total  machine  operating  time  that  the  system  will  achieve 
in  time  t  if  it  is  initialized  in  state  i.  Using  these  definitions,  it  follows  that 

Trj{t  +  dt)  =  Vj{t){l -'^Tjidt} +  '^v,{t)r,jdt  ji  =  l,2,...,/  (A.l) 

v,{dt  +  t)-{l-Yl  rtjdt){nl^Ut  +  ^>(0  *  =  1,2,...,/  ( A.2) 

J#' 
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Equations  A.l  and  A.2  can  be  simplified  by  defining  r„  =  —  and  substitut¬ 

ing  to  yield 


•Kj{t  -f  dt)  =  7rj(t){l  -f  Tjjdt}  +  ^  iri{t)r,jdt 

; 

1=1 


(A.3) 


Vi{t  +  dt)  =  (1  +  r„dt){no  +  ^«(0}  + 

/ 

=  rAl^dt  -f  u,(t)  -i-  no  -|-  ‘rtjVj{t)dt 

j=i 


(A.4) 


Simple  algebraic  manipulation  yields  the  expressions 


rj{t  +  dt)  -  irj{t)  ’ 

- -r. -  -  2^^i{i)r,j 


(A.5) 


V,(t  -h  dt)  Vi(t)  _  (i)  (i 


=  nj,‘'  -I-  n[,'V„dt  -t-  ^  rijVj{t) 

j=i 


(A.6) 


Taking  the  limit  dt  — ►  0  produces  the  constant-coefficient  difierential  equations 


dirj{t)  _  /  N 

j.  — 

•=! 

(A.7) 

^  =  n«+gr„»,(() 

i  =  l,2,...,/ 

(A.8) 

These  equations  relate  state  probabilities  to  transition  rates  and  govern  the  total 
expected  machine  operating  time  for  the  system  in  time  t. 

For  all  large  t  such  that  a  Markov  process  reaches  jsteady-state  operation, 
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it  is  known  that  dTj{t)/dt  =  0  for  each  state  j.  Thus,  Equations  A. 7  demand  that 

/ 

TT.T.j  =  0  j  =  1, 2, . . . ,  /  ( A.9) 

1=1 

where  each  tt,  is  a  steady-state  probability  for  all  large  t.  Further,  each  n,(t)  in  each 
Equation  A. 8  can  be  replaced  with  an  asymptotic  expression  E{no)t  +  r,,  where 
E(tio)  is  the  slope  of  the  asymptote  (rate  at  which  the  system  accumulates  machine 
operating  time)  and  v,  is  the  intercept.  Equations  A. 8  then  become 

/ 

E(no)  =  tiq  ^  -I-  53  rij{E{no)t  +  vj} 
j=i 

1  / 

=  no  ^  +  E{no)t  53  +  $3  (A.IO) 

j=i  i=i 

The  expression  “’'ij  is  zero  by  dehnition,  leaving  the  value  determination  equa¬ 
tions 


/ 

E{no)  =  njj’^  -I-  53 i  =  1,2,...,/  (A.ll) 

j=i 

This  set  of  J  equations  has  /  1  unknowns.  Each  v,  represents  the  relative  value 

of  beginning  in  a  certain  state,  so  any  one  of  them  can  be  set  to  an  arbitrary  value. 
For  consistency,  set  v/  =  0. 

Now  consider  two  manpower  assignment  policies  represented  by  decisions 
k'  and  k"  for  each  system  state.  Assume  that  the  policy  improvement  step  of 
the  algorithm  in  Section  2.2  has  produced  each  decision  k"  as  a  successor  to  each 
decision  k'.  It  is  therefore  known  that 


j=i  j=i 


(A.12) 
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or 


(A.13) 

j=i  i=i 

From  the  value  determination  Equations  A. 11, 

it  is  known  that 

E(no)"  =  no^  +  E’’o”^i 

j=i 

(A.14) 

E(no)'  =  4‘^  +  E4'^j 

i=i 

(A.15) 

Subtracting  Equations  A.  15  from  A.  14  yields 


EinoY'-Eino)'  = 

j=i  j=i 


J=1 


=  i.  +  E’-SK--;) 

j=i 


(A.16) 


Let  E(no)^  =  E(no)"  -  E(no)'  and  =  v"  -  v'.  Equations  A.16  then  become 

£!(no)^  =  ^. +  1=1,2,...,/  (A. 17) 

j=i 

Multiplying  each  Equation  A. 17  by  a  corresponding  steady-state  probability  tt"  and 
summing  over  all  i  yields 


1=1  1=1  j=i 

1=1  J=1  1=1 


(A.18) 


64 

However,  Equations  A. 9  require  that  =  0  for  each  state  j.  Thus, 

=  (A.19) 

1=1 

Since  all  tt"  >  0  and  all  Si  >  0,  it  follows  that  E{no)^  >  0.  Thus,  E{no)"  will  be 
greater  than  E(no)'  if  decisions  k"  can  be  found  such  that  Ylj=i 
for  any  state  i. 

It  remains  only  to  show  that  a  better  policy  cannot  exist  without  being 
found  at  some  time  by  the  policy  improvement  step.  Assume  that  E{no)"  >  E{nuy 
but  the  algorithm  has  converged  on  the  policy  represented  by  decisions  k'.  Then 
6,  <  0  for  all  i  when  all  Si  are  defined  by  Equations  A. 13.  Since  tt"  >  0  for  all  i, 
Equation  A.  19  demands  that  E(no)"  —  <  0.  However,  this  contradicts  the 

assumption  that  jE(no)"  >  E(no)'.  It  is  therefore  impossible  for  a  superior  policy 
to  exist  but  remain  undiscovered  at  the  completion  of  the  algorithm. 


Appendix  B 


Program  MAINTOP 


B.l  Description 

Maintenance  resource  optimization  using  the  sequential  linear  program¬ 
ming  approach  is  implemented  in  a  portable  computer  program  called  MAINTOP 
(MAINTenance  Optimization).  The  program  is  written  in  standard  Pascal  and  is 
listed  in  the  next  section  of  this  appendix.  The  listing  includes  file  statements  for 
both  mainframe  and  microcomputer  implementations.  Some  internal  documenta¬ 
tion  is  provided,  including  variable  definitions  and  brief  descriptions  for  each  sub¬ 
routine.  An  overview  of  the  relationship  between  the  various  program  components 
can  be  obtained  from  Figure  B.l. 

Program  MAINTOP  makes  extensive  use  of  recursive  calls  to  subroutines. 
This  approach  achieves  programming  efficiency  in  constructing  the  queuing  net¬ 
work,  building  the  system  state  space,  defining  alternative  resource  structures,  and 
identifying  nondominated  feasible  decisions  for  resource  assignments.  The  program 
also  employs  pointer  variables  and  linked  lists  to  preserve  computer  memory  and 
thus  increase  the  tractable  problem  size.  Maximum  use  of  computer  capability 
can  be  achieved  by  adjusting  the  constant  parameters  listed  at  the  beginning  of 
the  program.  During  program  execution,  network  reduction  will  be  performed  as 
necessary  to  limit  the  size  of  the  system  state  space. 

To  use  program  MAINTOP,  data  should  first  be  entered  into  a  text  file 
called  MAINTOP.DAT  in  accordance  with  the  format  displayed  in  Table  B.l.  The 


Figure  B.l:  Relationship  between  Components  of  Program  MAINTOP 
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general  format  is  paralleled  by  a  sample  data  file  for  the  “Mike’s  Flying  Club” 
problem  introduced  in  Chapter  2.  Precise  character  position  is  not  critical,  provided 
the  data  is  sequenced  properly  and  placed  in  separate  lines  for  each  maintenance 
task  and  resource.  For  each  task  m,  the  set  is  defined  as  the  set  of  all  tasks 
which  must  be  completed  before  task  m  can  be  initiated.  A  value  of  0  in  the 
Am  position  indicates  that  a  task  is  always  required  between  operational  activities. 
A  value  of  0  is  also  used  as  an  “end  of  list”  character  for  the  task  and  resource 
parameters.  Care  should  be  taken  to  avoid  trailing  blanks  on  any  line  of  data. 


Table  B.l:  Data  Format  for  Program  MAINTOP 


N  fio  C 


2  0.50  100.0 


m  flm  fflm  Am 

0 


1  1.00  1  0  2  3 

2  0.25  1  0.20 

3  0.50  2  0.25 

0 


y  Qy 
0 


1  10.00  1 
2  20.00  2 

3  25.00  3 

4  30.00  2  3 

5  33.00  1  2  3 

0 


Output  from  the  program  is  recorded  to  a  text  file  called  MAINTOP .  OUT.  The 
output  includes  an  echo  listing  of  input  parameters,  a  description  of  the  generated 
model,  and  a  presentation  of  performance  results.  Sample  output  for  the  “Mike’s 
Flying  Club”  problem  is  listed  in  the  final  section  of  this  appendix. 
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B.2  Program  Listing 


program  MAINTOP (input .output .Data, Out) ; 

■{*******m******mmmmm*mmmmmm*mm*mmm**m0*mmmm*mm*»mmm*mm*****m*mmmm*imm*mmmm*m»*m*y 


i*  *> 

{*  Determines  optimal  specialization  and  task  allocation  for  maintenance  *} 

{»  resources  using  a  queuing  network  formulation.  *} 

{*  *} 

{*  Authors:  Dennis  C.  Dietz  and  Matthew  Rosenshine  *} 

{*  Department  of  Industrial  and  Management  Systems  Engineering 

Pennsylvania  State  Oniversity,  University  Park,  PA  16802  *} 

*} 

{*  Language:  Steuidard  Pascal  —  program  is  suitable  for  mainframe  or  *} 

■{*  microcomputer  implementation.  *}- 

<•  *} 


const  MazTask=10; 
MaxRsrc=2S ; 
MaxSpec=100; 
MaxNode=250 ; 
MaxState^SOO; 
MaxDec=2000; 


{maximum  number 
{maximum  number 
{maximum  number 
{maximum  number 
{maximum  number 
{maximum  number 


of  task  t3^es} 

of  resource  types} 

of  specialization  strategies} 

of  network  nodes} 

of  system  states} 

of  decision  variables} 


type  TaskSet  =  set  of  l..MaxTask; 

RodeSet  -  set  of  O..MaxNode; 

RsrcSet  =  set  of  l..MaxRsrc; 

TaskType  =  record 
TaskRate:  real; 

ReqRsrc:  integer; 
FailRate:  real; 

Precede:  TaskSet 
end;  {TaskType} 

TaskArray  =  array  [1. 


{rate  of  task  completion} 

{number  of  maintenance  resources  required} 
{rate  of  random  failure} 

{set  of  tasks  this  task  must  precede} 


MaxTask]  of  TaskType; 
TaskIntArray  =  array  [1.. MaxTask]  of  integer; 
TaskTst Array  =  array  [1.. MaxTask]  of  boolean; 
RsrcType  =  record 


{unit  cost  of  resource} 

{set  of  tasks  this  resource  can  perform} 
{minimum  resource  availability} 

{maximum  resource  availability} 
{availability  increase  increment} 


Cost:  real; 

TaskQual:  TaskSet; 

MinRsrc:  integer; 

MaxRsrc :  integer ; 

StpRsrc :  integer 
end ;  {RsrcType} 

Rsrc Array  =  array  [l.. MaxRsrc]  of  R8rcT3rpe; 

RsrcIntArray  =  array  [1.. MaxRsrc]  of  integer; 

NodeType  =  record 

Prob:  real;  {routing  probability} 

RedAdj :  real;  {rate  adjustment  for  network  reduction} 

Pending:  TaskSet;  {set  of  pending  tasks} 


69 


Eligible:  TaskSet;  {set  of  eligible  tasks} 

NextModes:  NodeSet  {set  of  succeeding  nodes} 
end;  {NodeType} 

NodeArray  =  array  [O..HaxNode]  of  NodeType; 

NodeIntArray  =  array  [O..MaxNode]  of  integer; 

OccupyPointer  =  ‘OccupyType; 

OccupyType  =  record 

Node:  O..HaxNode;  {occupied  node  identification} 

Machs:  integer;  {nusiber  of  machines  at  node} 

NextOccupy:  OccupyPointer  {pointer  to  next  occupied  node} 
end;  {OccupyType} 

StateArray  =  array  [1 . .HaxState]  of  OccupyPointer; 
StateIntArray  =  array  [1 . .MaxState]  of  integer; 

StateRealArray  =  array  tl - -MaxState]  of  real; 

TransPointer  =  "TransType; 


{new  state  after  transition} 
{node  which  loses  a  machine} 
{node  which  gains  a  machine} 
{task  completed  at  transition} 
{pointer  to  next  transition} 


TransType  =  record 

NextState :  1 . . MaxState ; 

LossNode:  O..MaxNode; 

GainNode :  0 . -MaxNode ; 

TransTask :  1 . . MaxTask ; 

NextTrans:  TransPointer 
end;  {TransType} 

Trans Array  =  array  [1 . -MaxState]  of  TransPointer; 
DecMatrix  =  array  [0. .MaxState,0- -MaxDec]  of  real; 
ExPointer  =  "ExType; 

ExType  =  record 

ExDec:  1- -MaxDec; 

NextEx:  ExPointer 
end;  {ExType} 

AltPointer  =  "AltType; 

AltType  =  record 

Structure:  RsrcIntArray ; 

SpecCode :  1 . . HaxSpec ; 

Excluded:  ExPointer; 

CmtEx:  ExPointer; 

Accept:  boolean; 

GoodGuess:  real; 

Optimum:  real; 

NextAlt ,PrevAlt :  AltPointer 
end;  {AltType} 

Opt Array  =  array  [1- -HaxSpec]  of  AltPointer; 


{excluded  decision} 

{pointer  to  next  excluded  decision} 


{resource  structure} 

{specialization  strategy} 

{pointer  to  excluded  decisions} 

{pointer  to  current  excluded  decision} 
{indicator  for  acceptable  structure} 
{estimater  for  relative  performance} 
{expected  number  of  operating  machines} 
{pointer  to  next/previous  alternative} 


var  NufflMach:  integer; 
NumTask:  integer; 
NumRsrc:  integer; 
NumNode :  integer ; 
NufflState:  integer; 
NumDec :  integer ; 
NodeLimit:  integer; 
OpsRate:  real; 
CostLimit:  real; 


{number  of  machines  in  system} 

{nxunber  of  task  types} 

{number  of  resource  specialty  t3^es} 
{number  of  nodes  in  network} 

{number  of  system  states} 

{number  of  decision  variables} 

{number  of  nodes  in  reduced  network} 

{rate  of  operational  activity  completion} 
{resource  expenditure  limit} 
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Method:  char; 

TaskData:  TaskArray; 
LosNode:  Taskint Array ; 
RsrcData:  RsrcArray; 
NodeData:  NodeArray; 
StateData:  StateArray; 
TransData:  TransArray; 
AltData:  AltPointer; 

Dec:  DecMatrix; 

BasicDec:  StateIntArray ; 
Data, Out:  text; 


{solution  method} 

{array  of  task  data} 

{array  of  low  number  nodes  for  task  set  size} 
{array  of  resource  data} 

{array  of  node  data} 

{array  of  state  data} 

{array  of  transition  data} 

{pointer  to  list  of  structure  alternatives} 
{matrix  of  decision  coefficients} 

{array  of  basic  decisions} 

{input  and  output  text  files} 


procedure  ReadData(var  NumMach,NumTask,NumRsrc :  integer; 

var  OpsRate.CostLimit:  real; 
var  Method:  char; 
var  TaskData:  TaskArray; 
var  RsrcData:  RsrcArray); 


{ - > 

{  Reads  model  parameters  from  input  data  file  and  echos  to  output.  } 
{  } 
{  Calling  Routines:  program  MAINTOP  } 
{  Called  Routines:  none  } 
< - } 


var  I,J:  integer;  {local  counter  variables} 
begin 

writelnC’**  Reading  Data 

writeln(Out, INPUT  DATA  SUMMARY  ••••’); 
writeln(Out) ; 

readln(Data,N\unMach,OpsRate,CostLimit) ; 
vriteln(Out , ’Number  of  Machines:  ’ ,NumMach:8) ; 
writeln(Out , ’Operations  Rate:  ’ ,0psRate : 8 :4) ; 

writeln(Out , ’Expenditure  Limit:  ’ ,CostLimit : 8 : 2) ; 
write(0ut , ’Solution  Method:  ’); 

writeln; 

writeln(’ -  Solution  Method  Menu  - ’); 

writelnC ’  N.  Optimize  each  specialization  strategy  (NORMAL)’); 
writeln(  ’  I.  Evaluate  all  resource  structures  (INSIGHT)’); 

writeln(’  S.  Find  optimal  structure  only  (SPEED)’); 

writeln; 

Method: =’  ’; 
repeat 

write( ’Select  solution  method  — >  ’); 

readln(Method) ; 

writeln 

until  (Method  in  [’N’ , ’n’ , ’I’ , ’i’ , ’S’ , ’s’]) ; 
writeln; 
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case  Method  of 

>N’ , ’n’ :  writelnCOut, ’  MORMAL*): 

writeln(Out, >  IKSIGHT’): 

’S’, ’s’:  writeln(Out , ’  SPEED*) 
end;  -{case} 

WTiteln(Out) ; 
writeln(Out) ; 

NuniTask:=0; 

for  1:=1  to  HaxTask  do  TaskData [I] .Precede :=[] ; 

read(Data,I) ; 

while  (loO)  do  begin 

writelnCOut , ’Task  ’,1:2); 

NumTask :  =NuinTask-)-l ; 

with  TaskData [NumTask]  do  begin 

read(Data,TaskRate,ReqRsrc,FailRate) ; 
writeln(Out , ’Repair  Rate;  ’ ,TaskRate :8 ;4) ; 

writelnCOut , ’Resources  per  Task:  ’ ,ReqRsrc :8) ; 
write (Out,  ’Failure  Rate:  ’); 

if  (FailRate>0.0)  then  writeln(Out ,FailRate :8 :4) 
else  writelnCOut,’  — ’) 

end;  -{with} 

write(Out , ’Prerequisite  Tasks:  ’); 
while  not  eoln(Data)  do  begin 
read(Data, J) : 
write(Out , J :4) ; 

with  TaskData  [J]  do  Precede  :=Precede-*-tI] 
end;  {while} 
writeln(Out) ; 
writeln(Out) ; 
read(Data,I) 
end ;  {while} 
writeln(Out) ; 

NufflRsrc ; =0 ; 
read(Data,I) ; 
while  (loo)  do  begin 

writelnCOut , ’Resource  ’,1:2); 

NumRsrc : =NumRsrc+l ; 
with  RsrcDat a [NumRsrc]  do  begin 
read(Data,Co8t) ; 

writelnCOut , ’Unit  Cost;  ’,Cost:8;2); 

TaskQual : = [] ; 

write (Out, ’Qualified  Tasks:  ’); 
while  not  eoln(Data)  do  begin 
readCData, J) ; 
writeCOut , J :4) ; 

Task()ual :  =TaskQual-t-  [J] 
end;  {while} 
end;  {with} 
writeln(Out) ; 
writeln(Out) ; 
readCData, 1) 
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end;  {while} 
writeln(Out) ; 
end ;  {Re  adBat  a} 

^0jt*mmm**»***’t‘*m*mmm*mmm*m*mm*mmmmmmmm*t¥*mtt»*****m***»**m*m*mm**mmm*m***mm***’**y 

procedure  NakeNode(var  TaskSetSz ,NuniTask,NuinNode :  integer; 
var  OpsRate.RateSun:  real; 
var  TaskData:  TaskArray; 
v2Lr  LowNode:  Tasklnt Array ; 
var  NodeData;  NodeArray; 
var  TmpTaskSet .PrecedeSet :  TaskSet) ; 


^ - > 

{  Creates  a  network  node  and  computes  node  parameters .  } 
{  } 
{  Calling  Routines :  procedure  NextTask  } 
{  Called  Routines:  none  } 
^ - } 


var  I,J:  integer;  {local  counter  variables} 
begin 

NumNode ;  sNumNode-t-l ; 

if  (LowNode[Ta8kSetSz]=0)  then  LowNode[TaskSetSz] :=NuniNode; 
with  IlodeData[NumNode]  do  begin 
Pending : =TmpTa8kSet ; 

Eligible : =TmpTa8kSet-PrecedeSet ; 

Prob : sQpsRate/RateSum; 
for  I;=l  to  NusiTask  do 

if  (TaskData [I] .FailRate=0.0)  and  not  (I  in  Pending)  then  Prob:=0.0; 
if  (Prob>0.0)  then 

for  I:=l  to  NumNode-1  do 

if  (NodeDataEl] .Pending<=TmpTaskSet)  then 
Prob:=Prob-NodeData[I] .Prob; 

RedAdj :=1 .0; 

if  (Ta8kSetSz=l)  then  NextNodes : = [0] 
else  begin 

NextNodes :=[] ; 

for  I : =LowNode [TaskSetSz-l]  to  LowNodeCTaskSetSz]-!  do 
for  J;=l  to  NumTask  do 

if  (J  in  NodeData [NvmNode] .Eligible) 

and  (NodeDataEl] .Pending+ [J]=TmpTaskSet)  then 
NextNodes :  =NextN&des-i’  El] 
end;  {else} 
end;  {with} 
end;  {MakeNode} 

{***************4i**»**«**«*********»**************»4i***** ***••***«*****»»»**•**} 


procedure  HextTask(I,TaskSetSz:  integer; 
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var  HiFirstTask.TaskCnt:  integer; 

var  NumTask.NvunNode:  integer; 

var  OpsRate:  real; 

var  TaskData:  TaskArray; 

var  LonNode:  Taskint Array ; 

var  NodeData:  NodeArray; 

var  TfflpTaskSet:  TaskSet); 


{ - } 

{  Recursively  builds  a  pending  task  set  and  determines  if  a  node  with  the  } 
{  current  pending  task  set  should  be  created.  } 
{  } 
-[  Calling  Routines:  procedure  HakeNetuork,  procedure  NextTask  } 
{  Called  Routines:  procedure  HakeNode  } 
{ - } 


var  J:  integer; 

RateSum :  real ; 
PrecedeSet:  TaskSet; 
AluaysSet:  TaskSet; 


{local  counter  variable} 

{sum  of  eligible  task  completion  rates} 
{tasks  which  must  precede  all  pending  tasks} 
{pending  tasks  which  are  always  required} 


begin 

if  (NumNode<MaxNode)  then  begin 
TmpTaskSet : =TmpTa8kSet+ [I] ; 

TaskCnt : =Ta8kCnt+l ; 
if  (TaskCnt^TaskSetSz)  then  begin 
PrecedeSet :  =  [] ; 

AlwaysSet:-[] ; 

RateSum : =0p8Rate ; 
for  J:=l  to  NumTask  do 

if  (J  in  TmpTaskSet)  then 

PrecedeSet : =PrecedeSet+TaskData [J] . Precede 
else  RateSum  :=RateSun-t’Ta8kData[J]  .FailRate; 
for  J:=l  to  NvunTask  do 

if  (J  in  PrecedeSet)  and  (TaskData[J] .FailRate=0.0)  then 
AlwaysSet :  =AlwaysSet-'’  [J]  ; 
if  ( AlwaysSet <=TmpTa8kSet)  then 

HakeNode ( TaskSetSz , NumTask , NumNode , OpsRate , RateSum .TaskData , 
LowNode , NodeData , TmpTaskSet , PrecedeSet ) 
end  else  begin 

HiFirstTask : =HiFir8tTa8k+l ; 
for  J:=I+1  to  HiFirstTask  do 

NextTask( J .TaskSetSz, HiFirstTask, TaskCnt .NumTask, NumNode , 
OpsRate .TaskData , LowNode .NodeData , TmpTaskSet ) ; 

HiFirstTask : =HiFir8tTa8k-l 
end;  {else} 

TmpTaskSet :=TmpTa8kSet- [I]  ; 

TaskCnt :=Ta8kCnt-l 
end; 

end;  {NextTask} 


74 


■^mm*mm***m^‘m**m****m*m****0mmmmm*m****mmmm*»*********m***»mmtr*****m*m*mm0*<¥mm*^y 

procedure  MakeMetwork(var  NumTask,NumNode:  integer; 

var  OpsRate:  real; 
var  TaskData:  TaskArray; 
var  LowNode:  TaskIntArray ; 
var  NodeData:  N ode Array ) ; 


{ - 

{  Creates  a  queuing  network  from  model  parameters. 

{ 

{  Calling  Routines:  program  MAINTOP 

{  Called  Routines:  procedure  NextTask 

{ - 


■} 

} 

} 

} 

} 

■> 


var  I:  integer; 

TmpTaskSet :  TaskSet ; 
HiFirstTask:  integer; 
TaskSetSz:  integer; 
TaskCnt :  integer ; 


-[local  counter  vao-iable} 

■[temporary  set  of  pending  tasks} 
{highest  numbered  first  task  in  set} 
{number  of  tasks  in  set} 

{task  counter} 


begin 

Hriteln(’**  Building  Network 

for  I:=l  to  NumTask  do  LowNode [1] :=0; 

NumNode : =0 ; 

TmpTaskSet : = [] ; 

TaskCnt :=0; 

for  TaskSetSz :=1  to  NumTask  do  begin 
HiFirstTask :  ^NumTask-TaskSetSz-i-l ; 
for  I:=l  to  HiFirstTask  do 

NextTask( I , TaskSetSz .HiFirstTask , TaskCnt , NumTask , NumNode , 
OpsRate , TaskData , LowN  ode , NodeData , TmpT askSet ) 
end;  {for} 
end;  {MakeNetwork} 


function  NodeCapCvar  NumNach, NumNode :  integer):  integer; 


{ - } 

{  Returns  the  maximum  number  of  nodes  permitted  by  the  maximum  system  } 
{  state  size  (HaxState) .  } 
{  > 
{  Calling  Routines:  program  MAINTOP  } 
{  Called  Routines:  none  } 
{ - } 


var  I:  integer;  {local  counter  variable} 

StSpaceSz:  integer;  {state  space  size} 
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begin 

StSpaceSz :  =NumNach-t-l ; 

I:=l: 

vhile  (I<=NumNode)  and  (StSpaceSz<HaxState)  do  begin 
I: =1+1; 

StSpaceSz ;=StSpaceSz*(NuiDMach+I)  div  I 
end;  ■[ while} 

NodeCap:=I-l 
end;  {NodeCap} 


procedure  ReduceNetworkCvar  NumTask.NumNode.NodeLinit:  integer; 

var  TaskData:  TaskArray; 
var  LowNode:  Taskint Array ; 
var  NodeData:  NodeArray) ; 


{ - > 

{  Reduces  network  size  and  adjusts  node  parameters.  } 

i  > 

{  Calling  Routines :  program  MAINTOP  } 

-[  Called  Routines:  none  } 

{ - } 


var  integer; 

TaskSetSz:  integer; 
Ratel:  real; 

RateJ:  real; 

ProblJ :  real; 
NewProbJ:  real; 
NewRateJ :  real; 


{local  counter  variables} 

{number  of  tasks  in  set} 

{departure  rate  from  losing  node} 

{departure  rate  from  gaining  node} 

{transition  probability  from  losing  to  gaining  node} 
{adjusted  routing  probability} 

{adjusted  departure  rate} 


begin 

writelnC ’ •*  Reducing  Network 
for  I:=NumNode  downto  NodeLimit+1  do  begin 
TaskSetSz : =0 ; 

Ratel; =0.0; 

for  K:=l  to  NumTask  do  begin 

if  (K  in  NodeData [1] .Pending)  then  TaskSetSz : =Ta8kSetSz+l ; 
if  (K  in  NodeOataCi] .Eligible)  then 
Ratel :=RateI+TaskData[K] .TaskRate 
end;  {for} 

for  J:=LowNodeCTaskSetSz-l]  to  LowNode [TaskSetSz] -1  do 
if  (J  in  NodeData [I] .NextNodes)  then  begin 
RateJ: =0.0; 

for  K:=l  to  NumTask  do  begin 

if  (K  in  NodeDataCJ] .Eligible)  then 
RateJ :=RateJ+Ta8kData[K] .TaskRate; 
if  (NodeDataCJ] .Eligible+[K]=NodeData[I] .Eligible)  then 
ProblJ :=Ta8kData[K] .TaskRate/Ratel 
end;  {for} 
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NewProbJ :=NodeData[J] .Prob+NodeData[I] .Prob*ProbIJ ; 
NenRateJ:=NeBProbJ/(NodeData[J] .Prob/(RateJ*NodeData[J] .RedAdj) 
+NodeData[I] .Prob*ProblJ* (1.0/ (Rate J*NodeData[J] .RedAdj) 

+1 .0/(RateI*NodeData[I] .RedAdj) ) ) ; 

NodeData[J] .Prob;=NesProbJ ; 

NodeData[J]  .RedAdj  :=NeHRateJ/(RateJ'*NodeData[Jj  .RedAdj) 
end;  {if} 
end;  {for} 

end;  {ReduceNetvork} 

{«»«»«»»»»**«*4i*#**4i*W»«****««4i4>**«»**i»»»**i«**»«4i**»**«*»****i|i«*»»***4i*««»»«*<*#} 

procedure  ListNet«ork(var  NumTask.NumNode.NodeLimit:  integer; 

var  NodeData:  NodeArray); 


.[ - } 

{  Writes  network  parameters  to  output  file.  } 
{  } 
{  Calling  Routines:  program  MAINTOP  } 
{  Called  Routines:  none  } 
.[ - > 


var  I,J:  integer;  {local  counter  variables} 
begin 

Hriteln(0ut, ’••••  MODEL  SUMMARY  **••’); 
writeln(Out) ; 

writeln(0ut , ’Number  of  maintenance  nodes  in  full  network:  ’ ,NumNode:4) ; 
if  (NodeLimit<NumNode)  then 

writeln(0ut , ’Number  of  maintenance  nodes  in  reduced  network:  ’, 
NodeLimit :4) ; 

if  ( NodeLimit cNumTask)  then 

writeln(0ut , ’WARNING :  Fewer  allowable  nodes  than  tasks.’); 
writeln(Out) ; 

for  I:=l  to  NodeLimit  do  begin 
writeln(0ut , ’Node  ’,1:4); 
with  NodeData [I]  do  begin 

writeln(0ut , ’Routing  Probability:  ’,Prob:5:4); 

if  (NodeLimit<NumNode)  then 

writeln(0ut , ’Reduction  Adjustment:  ’, RedAdj :5:4) ; 
write(0ut, ’Pending  Tasks:  ’); 
for  J:=l  to  NumTask  do 

if  (J  in  Pending)  then  write ( Out ,J:4) ; 
writeln(Out) ; 

write ( Out , ’Eligible  Tasks :  ’); 
for  J:=l  to  NumTask  do 

if  (J  in  Eligible)  then  write(0ut,J:4) ; 
writeln(Out) 
end;  {with} 
writeln(Out) 
end;  {for} 
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end;  {ListNetwork} 

■{»Lmm**‘mmm**»*m*m*»m*mmm**mm*****l***mm*mmm«**mm**mmm**mmm***m^‘*m*m**^^0^*^^**m*•mm*y 

procedure  MadceState(v2u:  NodeLimit .NumState:  integer; 

var  StateData:  StateArray; 
var  TmpState:  ModeIntArray) ; 


^ - } 

{  Creates  a  system  state.  The  machine  distribution  is  stored  in  a  linked  } 
■{  list  of  occupied  nodes .  } 
{  } 
{  Calling  Routines:  procedure  NextTrans  } 
{  Called  Routines;  none  } 
{ - } 


var  I:  integer;  {local  counter  variable} 

CmtOccupy:  OccupyPointer ;  {pointer  to  current  occupied  node} 

begin 

NumState :  =NumState-t-l ; 
for  I:-0  to  NodeLimit  do 

if  (TmpState [I] >0)  then  begin 

if  (StateData [NtuRState}=nil)  then  begin 
nea(StateOata[NximStateJ ) ; 

CmtOccupy ;  =StateOata  [NumState] 
end  else  begin 

new  (CmtOccupy*  .NextOccupy) ; 

CmtOccupy :  =CmtOccupy*  .NextOccupy 
end;  {else} 

with  CmtOccupy*  do  begin 
Node:=I; 

Hachs ; =7mpState [I] ; 

NextOccupy :  mil 
end;  {with} 
end;  {if} 
end;  {NakeState} 

{«»*******«*««*«««««***«*i»«***««**«4t***«***********«»****»***********» »**»*•*•«} 

procedure  NextNoded:  integer; 

var  NodeCnt, NodeLimit, NumState:  integer; 
var  StateData:  StateArray; 
var  TmpState:  NodelntArray) ; 


< - ] 

{  Recursively  distributes  machines  to  network  nodes  and  identifies  system  ] 
{  states.  ] 
{  ] 
{  Calling  Routines:  procedure  MakeStateSpace,  procedure  NextNode  ] 
{  Called  Routines:  procedure  NakeState  ] 
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{ - } 

var  J:  integer;  {local  counter  variable} 
begin 

NodeCnt :  =NodeCnt-i-l ; 
if  ( NodeCnt =NodeLimit)  then  begin 
TmpState [NodeCnt] : =I ; 

HakeStateCNodeLinit .NimState .StateDat a, TmpState) 
end  else 

for  J:=0  to  I  do  begin 

TmpState [NodeCnt] : =I-J ; 

Next  N  ode ( J , N  odeCnt , N  odeLimit , NumState , St at eData , TmpState ) 
end;  {for} 

NodeCnt : =NodeCnt-l 
end;  {NextNode} 

procedure  MakeStateSpaceCvar  NumMach.NodeLimit .NumState;  integer; 

var  StateData:  StateArray) ; 


i - > 

{  Creates  the  system  state  space.  } 
{  } 
{  Calling  Routines:  program  MAINTOP  } 
{  Called  Routines:  procedure  NextNode  } 
i - } 


var  I:  integer;  {local  coiinter  variable} 

NodeCnt;  integer;  {node  counter} 

TmpState:  Nodelnt Array ;  {temporary  state  data} 

begin 

vriteln(’**  Defining  System  States  _ ’); 

for  I:=l  to  HaxState  do  StateData[I] :=nil; 

NumState :=0; 

NodeCnt :=0; 

for  I:=0  to  NumNach  do  begin 

TmpState [NodeCnt] : =NumMach-I ; 

NextNoded .NodeCnt .NodeLimit .NumState, StateData, TmpState) 
end;  {for} 

uritelnCOut , ’Number  of  system  states:  ’ .NumState :6) ; 
sritelnCOut) 
end;  {HakeStateSpace} 

{«««*«««*«*«««4t«**«**4i*«»***«**«***«*******«4i«*«**i»*i»*****«»««**»****«»»*«*»»**}' 

procedure  MakeTrans ( I , J ;  integer; 

var  Loser .Gainer .lumTask:  integer; 
var  TaskData:  TaskArray; 


var  NodeData:  NodeArray; 
var  TransData:  TransArray; 
var  CmtTrans:  TransPointer) ; 
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< - 

-{  Creates  a  data  record  for  transition  parameters. 

<. 

Calling  Routines:  procedure  MaReTransSpace 
■{  Called  Routines:  none 

{ - 


•} 

} 

} 

} 

> 

•} 


var  K:  integer;  {local  counter  variable} 

TmpTaskSet:  TaskSet;  {temporary  set  of  pending  tasks} 

begin 

if  (TransData [I] =nil)  then  begin 
ne«(TransData[I] ) ; 

CmtTrans :  =TransData  [I] 
end  else  begin 

new  ( CmtTrans  * .  NextTrans ) ; 

CmtTrans :  =CmtTrans *  .NextTrans 
end;  {else} 

with  CmtTrans*  do  begin 
NextState:=J ; 
if  (LosersO)  then  begin 
LossNode:=0; 

GainNode : =Gainer ; 

TransTask:=0 
end  else  begin 
LossNode : =Lo8er ; 

GainNode : =Gainer ; 

if  (Gainer=0)  then  TmpTaskSet :=[] 

else  TmpTaskSet ; =NodeData [Gainer] . Pending ; 
for  K:=l  to  NumTask  do 

if  (  [K]-»^TmpTa8kSet=NodeData [Loser]  .Pending)  then  TransTask:=K 
end;  {else} 

NextTrans ;  mil 
end;  {with} 
end;  {HakeTrans} 


.{«****»««««[««*«*»4>««*»«*«*»**«*«**«*«*«»*****«*******»«******«*»»«**«*«*4<*«^*i»**} 

procedure  MakeTransSpace(var  NumTask, NumState.NodeLimit:  integer; 

var  NodeData:  NodeArray; 
var  StateData:  StateArray; 
var  TransData:  TransArray); 


^ - } 

{  Creates  a  table  of  transition  data.  For  each  state,  transition  data  is  } 
{  stored  in  a  linked  list  of  records  to  save  memory.  } 
{  > 


{  Calling  Routines:  program  MAINTOP 

{  Called  Routines:  none 

^ - 
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} 

} 

•} 


var  I,J,K:  integer; 

ChgMach:  integer; 

Loser , Gainer:  integer; 

TmpStatel ,TmpStateJ :  NodeIntArray ; 

CmtI.CmtJ:  OccupyPointer ; 

CmtTrans:  TransPointer ; 

OneStep:  boolean; 

begin 

writeln(’**  Defining  Transitions 
for  I:=l  to  NufflState  do  begin 
TransData[I] :=nil; 

for  K:=0  to  Nodelimit  do  TmpStatel [K] :=0; 

CmtI:=StateDataCl]  ; 

Bhile  (CmtIOnil)  do  begin 

TmpStatel  [Cmtl*  .Node]  :  =CmtI*  .Machs ; 

CmtI:=CmtI*  .NextOccupy 
end;  {while} 

for  J:  =  l  to  NumState  do  if  (loj)  then  begin 
for  K:=0  to  NodeLimit  do  TmpStatel [K] :=0; 

CmtJ:=StateData[J]  ; 
while  (Cmtlonil)  do  begin 

TmpStatel  [Cmtl*  .Node]  ;sCmtl*  .Machs ; 

Cmtl:=Cmtl*  .NextOccupy 
end;  {while} 

OneStep : =true ; 

Loser :=-l; 

Gainer :=-l; 

K:=0; 

while  (OneStep=true)  and  (K<=NodeLimit)  do  begin 
ChgMach : =TmpStatel [K] -TmpStatel [K]  ; 
if  (abs (ChgMach) >1)  then  OneStep : =falBe 
else  case  ChgMach  of 

-1:  if  (Loser>=0)  then  OneStep :=false  else  Loser :=K; 

1:  if  (Gainer>=0)  then  OneStep :=false  else  Gainer :=K 
end;  {case} 

K:=K+1 

end ;  {while} 
if  (OneStep=true)  and 

((Loser=0)  or  (Gainer  in  NodeData[Loser] .NextNodes))  then 

MakeTr ans ( I , 1 , Loser , Gainer , NumTask , TaskData , NodeData , TransData , 
CmtTrans ) 
end;  {for} 
end;  {for} 

end;  {MakeTransSpace} 


{local  counter  variables} 

{change  in  machines  at  a  node} 
{losing  and  gaining  nodes} 
{temporary  state  data} 

{pointers  to  current  state  data} 
{pointer  to  current  trusition} 
{indicator  for  one-step  behavior} 
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procedure  MakeAlt(var  SpecCnt:  integer; 

var  AltData.CmtAlt :  AltPointer; 
var  TmpStructure :  Rsrcint Array ) ; 


{ - y 

■{  Creates  a  maintenance  resource  structure  alternative.  } 
{  } 
-C  Calling  Routines :  procedure  NextAvlRsrc  } 
{  Called  Routines:  none  } 
{ - } 


begin 

if  (AltData=nil)  then  begin 
new(AltData) ; 

AltData* . PrevAlt : =nil ; 

Cmt  Alt :  =AltData 
end  else  begin 

neH(CmtAlt*  .HextAlt) ; 

Cmt  Alt  *  .  Next  Alt  *  .  PrevAlt :  =Crat  Alt ; 

Cmt  Alt :  =Cmt  Alt  *  .  Next  Alt 
end;  {else} 

with  Cmt  Alt*  do  begin 

Structure : =TmpStructure ; 

SpecCode : =SpecCnt ; 

Excluded :=nil; 

CmtEx:=nil; 

NextAlt :=nil 
end ;  {sith} 
end;  {NakeAlt} 

{***««**i><)>»*<**»***«»«»4>»****4i«»**4>**»««i*****»***»««*»*»*4>*«>************«***»***} 

procedure  NextAvlRsrc (I:  integer; 

CostRem:  real; 

var  SpecCnt, NuiiiRsrc:  integer; 
var  RsrcData:  RsrcArray; 
var  AltData,CmtAlt :  AltPointer; 
var  TmpStmcture:  Rsrc Int Array  )  ; 


{ - } 

{  Recursively  adds  available  maintenance  resources  to  create  stmcture  } 
{  alternatives.  } 
{  } 
{  Calling  Routines:  procedure  AddAlts,  procedure  NextAvlRsrc  } 
{  Called  Routines:  procedure  NakeAlt  } 
^ - } 


var  J:  integer;  {local  counter  variable} 

var  Accept:  boolean;  {indicator  for  acceptable  alternative} 
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begin 

while  (I<=NuinRsrc)  and  (TmpStructure[I]=0)  do  I:=I+1; 
if  (I>IiuinRsrc)  then  begin 

if  (CostReffl>=0 . 0)  then  begin 
Accept : =true ; 
for  J:=l  to  NumRsrc  do 
with  RsrcDataCJ]  do 

if  (TmpStructure[J]>0)  and  (MaxRsrc-TmpStmctuTe[J]>=StpRsrc) 
and  (Co8t*StpRsrc<=CostRem)  then  Accept :=false; 
if  Accept  then  HaheAlt(SpecCnt,AltData,CmtAlt,TmpStructiire) 
end;  -(if} 
end  else 

with  RsrcDataCi]  do  begin 
J : =MinRsrc ; 

while  (J<=NaxRsrc)  do  begin 

if  (Co8tReni>=0.0)  then  begin 
TmpStructure [I] : = J ; 

NextAvlRsrcCI-t-l  ,CostRem,SpecCnt  .NumRsrc  .RsrcData, 

AltData , Cmt Alt  .TmpStructure ) 
end;  -(if} 

J : =J+StpR8rc ; 

CostRem : =CostRem-StpR8rc*Cost 
end;  {while} 
end;  {with} 
end;  {NextAvlRsrc} 

{«•*•*•••••«•*««««***•««••*»«**«•«**•***•«»•••*•«>*•* **«•»»*»•»»••*•••*• »*•«***•} 

procedure  AddAltsCvar  SpecCnt.NumNach.NumTask.NodeLimit.NiusRsrc:  integer; 
var  CostLimit:  real; 
var  TaskData:  TaskArray; 
var  NodeData:  NodeArray; 
var  RsrcData:  RsrcArray; 
var  AltData.  Cmt  Alt:  AltPointer; 
var  TmpStructure:  RsrcIntArray) ; 


< - > 

{  Adds  resource  structure  alternatives  to  a  linked  list.  } 
{  } 
{  Calling  Routines:  procedure  NextRsrc  } 
{  Called  Routines:  procedure  NextAvlRsrc  } 
^ - } 


var  I.J.K:  integer;  {local  counter  variables} 

NodeReq:  integer;  {number  of  resources  required  at  node} 
CostRem:  real;  {remaining  available  expenditure} 

begin 

CostRem : =CostLimit ; 

SpecCnt :  =SpecCnt-i-l ; 
write(0ut.SpecCnt:8,’  {  ’); 
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for  I:=l  to  NumRsrc  do  begin 

if  (TmpStmcture[I]>0)  then  begin 
vrite(0ut ,I : 1 , ’  ’): 
vith  RsrcData[I]  do  begin 
MinRsrc:=0; 

MaxRsrc : =0 ; 

StpR8rc:=9999; 

for  J:=l  to  NunTask  do 

if  (J  in  TaskQual)  then  begin 

if  (TaskData[J] .ReqRsrc>MinRBrc)  then 
HinR8rc:=TaskDataCJ] .ReqRBrc; 
if  (Ta8kData[J] .ReqR8rc<StpR8rc)  then 
StpR8rc:=Ta8kDataCJ] .ReqRsrc 
end;  {if} 

for  J:=l  to  NodeLimit  do  begin 
NodeReq:=0; 

for  K:=l  to  NunTask  do  begin 

if  (K  in  TaskQual)  and  (K  in  NodeData[J] .Eligible)  then 
NodeReq :  =NodeReq-t’TaskData  [K]  .  ReqRsrc *NunMach ; 
if  (NodeReq>MaxRsrc)  then  MaxRsrc : =NodeReq 
end;  {for} 
end;  {for} 

CostRem : =Co8tRem-MinRsrc*Co8t 
end;  {with} 
end;  {if} 
end;  {for} 
writeln(0ut, ’}’) : 

I:=l; 

Next AvlRsrc ( I , CostRem , SpecCnt , NumRsrc , RsrcData , AltData , Cmt Alt .TmpStructure ) 
end;  {AddAlts} 

{a******************************************************* ****••*****•*****•****} 

procedure  NextRsrcd:  integer; 

TaskCover:  TaskTst Array ; 

var  SpecCnt ,NunNach, NunTask, NodeLimit .NumRsrc:  integer; 

var  CostLimit;  real; 

var  TaskData:  TaskArray; 

var  NodeData:  NodeArray; 

var  RsrcData:  RsrcArray; 

var  AltData,  Cmt  Alt:  AltPointer; 

var  TmpStructure :  RsrcIntArray) ; 


{ - } 

{  Recursively  determines  acceptable  resource  specialization  strategies.  } 
{  Adds  alternative  resource  structures  whenever  an  acceptable  strategy  } 
{  is  found.  } 
f  } 
{  C.-lling  Routines:  procedure  NakeAltSpace ,  procedure  NextRsrc  } 
{  Called  Routines:  procedure  AddAlts  } 
{ - } 
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var  J:  integer;  {local  counter  variable} 

Cont:  boolean;  {indicator  for  continued  addition  of  resource  types} 
Accept:  boolean;  {indicator  for  acceptable  strategy} 

begin 

TmpStructure [I] : =1 ; 

Cont : =true ; 

for  J:=l  to  NumTask  do 

if  (J  in  RsrcDataEl] .TaskQual)  then  begin 

if  TaskCover[J]  then  Cont:=fal8e  else  TaskCover[J] :=true 
end;  {if} 

if  Cont  then  begin 
Accept : =true ; 
for  J:=l  to  NumTask  do 

if  not  TaskCoverEJ]  then  Accept :=false ; 
if  Accept  then 

AddAlt s ( SpecCnt , NumNach , NumTask , N  odeLimi t , NumRsr c , CostLimit , 
TaskData ,  N  odeData ,  Rsr  cData ,  AltDat  a ,  Cmt  Alt ,  TmpStructure ) 
else  for  J:=I+1  to  NumRsrc  do 

NextRsrc ( J , TaskCover , SpecCnt ,NumHach , NumTask , NodeLimit , NumRsrc , 
CostLimit , TaskData, NodeData,RsrcData,AltData,CmtAlt , 
TmpStructure) ; 

end;  {if} 

TmpStructtire  [I]  :  =0 ; 
end;  {NextRsrc} 

procedure  NakeAltSpaceCvar  NumNach, NumTask, NodeLimit, NumRsrc:  integer; 

var  CostLimit:  real; 
var  TaskData:  TaskArray; 
var  NodeData:  NodeArray; 
var  RsrcData:  Rsrc Array; 
var  AltData:  AltPointer) ; 


^ - 

{  Creates  a  linked  list  of  resoiirce  structure  alternatives. 

{ 

{  Calling  Routines:  program  MAINTOP 

{  Called  Routines:  procedure  NextRsrc 

{ - 

var  I , J :  integer ; 

SpecCnt :  integer ; 

TaskCover:  TaskTstArray ; 

TmpStructure :  Rsrcint Array ; 

Cmt  Alt:  AltPointer; 

begin 

vritelnC**  Building  Resource  Structures 


{local  counter  variables} 

{counter  for  specialization  strategy} 
{indicator  for  covered  tasks} 
{temporary  resource  struct\ire} 
{pointer  to  current  alternative} 


■} 

} 

} 

} 

} 

■} 
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AltData:=nil; 

SpecCnt : =0 ; 

lor  I;=l  to  NumRsrc  do  TmpStructure[I] ;=0; 
writelii(0ut , ’Specialization  Resources ’ ) ; 
writeln(0ut , ’  Strategy  Employed’); 

for  I:=l  to  NumRsrc  do  begin 

for  J:=l  to  NumTask  do  TaskCover [J] :=f alse ; 

NextRsrc ( I .TaskCover .SpecCnt ,NumMach,NumTask,NodeLimit .NumRsrc . 

CostLimit . TaskData . NodeData .RsrcDat a . AltData . Cmt Alt .TmpStructure ) 
end;  {for} 
writeln(Out) ; 

Hriteln(Out) 
end;  -{NakeAltSpace} 


procedure  MakeDec(var  I.NumDec:  integer; 

var  TmpState;  NodeIntArray ; 
var  AltData:  AltPointer; 
var  TmpDec:  StateRealArray ; 
var  Dec ;  DecMatrix ; 
var  BasicDec:  StateIntArray) ; 


{ - - 

{  Creates  a  nondominated  feasible  resource  assignment  decision. 

{ 

{  Calling  Routines:  procedure  NextTrans 

{  Called  Routines :  none 

{ . . . . 

var  J:  integer;  {local  counter  variable} 

CrntAlt:  AltPointer;  {pointer  to  current  alternative} 

begin 

NumDec :  =N\imDec+l ; 

if  (BasicDec [I] =0)  then  BasicDec [I] :=NumDec; 

Dec [0 .NumDec] : =TmpState [0] ; 

Dec [I .NumDec] :=0.0; 

for  J;  =  l  to  NumState  do  if  (loj)  then  begin 
Dec [J .NumDec] : =-TmpDec [J] ; 

Dec [I .NumDec] :=Dec [I. NumDec] +TmpDec[J] 
end;  {for} 

Cmt  Alt :  =  AltData ; 

while  (CmtAltOnil)  do  begin 

if  not  CmtAlt' .Accept  then  begin 

if  (CmtAlt*  .Excluded=nil)  then  begin 
new  (CmtAlt*  .Excluded); 

CmtAlt*  . CmtEx  :=CmtAlt*  .Excluded 
end  else  begin 

new(CmtAlt*  .CmtEx*  .NextEx) ; 

CmtAlt*  .CmtEx :=Crat Alt*  .CmtEx*  .NextEx 


•} 

} 

} 

} 

} 

•} 
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end;  {else} 

with  CrntAlt*  .CmtEx"  do  begin 
ExDec :  =Nu]nDec ; 

MextEx:=nil 
end ;  {with} 
end;  {if} 

CmtAlt  :=CrntAlt*  .NextAlt 
end;  {while} 
end;  {MakeDec} 

■{m*’***i^mm*itiimmmmm*mm*mmmm*m**mmm*<*m»**m**m»m*m***m*******<^»****m**’¥*mm***<Hf****my 

procedure  NextTrans(I:  integer; 

var  CmtTrans:  TransPointer ; 

var  TmpState:  NodelntArray ; 

var  AssnTask.HaxAssnTask:  Tasklnt Array ; 

var  TmpDec:  StateRealArray ; 

var  NufflTask.NumRsrc.NumState:  integer; 

var  OpsRate:  real; 

var  TaskData:  TaskArray; 

var  RsrcData:  RsrcArray; 

var  NodeData:  NodeArray; 

var  AltData:  AltPointer; 

var  Dec:  DecMatrix; 

var  BasicDec:  Stateint Array ) ; 


< - } 

{  Recursively  considers  possible  transitions  to  build  an  assignment  } 
{  decision.  } 
{  } 
{  Calling  Routines:  procedure  HakeDecSpace ,  procedure  NextTrans  } 
{  Called  Routines:  procedure  HakeDec  } 
^ - } 


var  J,K,L:  integer; 
Accept:  boolean; 
CmtAlt:  AltPointer; 
QnassnRsrc:  integer; 


{local  counter  variables} 
{indicator  for  acceptable  decision} 
{pointer  to  current  alternative} 
{resources  not  assigned  to  tasks} 


begin 

if  (CmtTrans=nil)  then  begin 
Accept :=false; 

CmtAlt :  =AltData; 
while  (CmtAltonil)  do  begin 
CmtAlt  *  .  Accept :  =tme ; 
for  J:=l  to  NumRsrc  do 

if  (CmtAlf  .Stmcture[J]>0)  then  begin 
DnassnRsrc:=CmtAlt*.Stmcture[J]  ; 
for  K:=l  to  lumTask  do 

if  (K  in  RsrcDataCJ] .TaskQual)  then 

OnassnRsrc : sDnassnRsrc- AssnTask [K] ^TaskData [K] . ReqRsrc ; 
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if  (UnassiLB.src<0)  then  CmtAlt*  .Accept  :=false 
else  for  K:=l  to  NumTask  do 

if  (K  in  RsrcData[J] .TaskQual)  and 
( AssnTask [K] <MaxA8snTask[K] )  and 

(UnassnRsrO^TaskDataEK]  .ReqRsrc)  then 
CmtAlt*  .Accept:=f  alse 

end;  {if} 

if  CmtAlt*  .  Accept  then  Accept :  =true ; 

CmtAlt :  =CmtAlt*  .NextAlt 
end;  {while} 
if  Accept  then 

MakeDec ( I , NumDec ,TnipState , AltData.TmpDec ,Dec .BasicDec) ; 
end  else  begin 

J :  =CmtTrans  *  .  NextState ; 
if  (CmtTrans*  .Lo8sNode=0)  then  begin 

TmpDec [J] ;=TmpState [0]*0psRate*NodeData[CmtTran8" .GainNode] .Prob; 
NextTrans (I , CmtTrans"  .NextTrans  ,TinpState, AssnTaek.MaxAssnTaak , 
TmpDec , NumTask , NumRsrc , NumState , OpsRate , TaskData ,RsrcDat a , 
NodeData,AltData,Dec,BasicDec) ; 
end  else  begin 

L :  =CmtTrans  *  .  TransTask ; 

for  K:=TmpState [CmtTrans*. LossNode]  downto  0  do  begin 
TmpDec [J] :=K*TaskData[L] .TaskRate 

♦NodeData [CmtTrans*  .LossNode]  .RedAdj  ; 

AssnTask  [L]  :=AssnTask[L]-«^K: 

NextTrans ( I , CmtTrans  * . NextTrans , TmpState .AssnTask .Max AssnTask . 
TmpDec . NumTask . NiuoRsrc . NumState . OpsRate . TaskData .RsrcDat a . 
NodeData.AltData.Dec.BasicDec) ; 

AssnTask [L] : -AssnTask [L]-K 
end;  {for} 
end;  {else} 
end;  {else} 
end;  {NextTrans} 

{****4i4i*^**4(*4ii*i*i^****i^****«***********4t^****4>*>l>********>^********4t***«*4t*4E*«***i^} 

procedure  MakeDecSpace( var  NumTask. NumRsrc. NodeLimit. NumState, NumDec;  integer; 

var  OpsRate:  real; 
var  TaskData:  TaskArray; 
var  NodeData:  NodeArray; 
var  RsrcData:  RsrcArray; 
var  StateData:  StateArray; 
var  AltData:  AltPointer; 
var  Dec:  DecNatrix; 
var  BasicDec:  StateIntArray) ; 


.[ - } 

{  Determines  nondominated  feasible  decisions  for  each  system  state  and  } 
{  resource  stmcture .  } 
{  } 
{  Calling  Routines :  procedure  Optimize  } 
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■[  Called  Routines:  procedure  NextTrans  } 

^ - > 


var  I,J;  integer;  {local  counter  variables} 

TmpState:  NodeIntArray ;  {temporary  state  data} 

TmpDec :  StateRealArray ;  {temporary  rate  data} 

AssnTask:  Taskint Array ;  {number  of  tasks  assigned  resources} 

MaxAssnTask:  Taskint Array ;  {maximum  nvimber  of  assignable  tasks} 

CmtOccupy:  OccupyPointer ;  {pointer  to  current  occupied  node} 

CmtTrans:  TransPointer ;  {pointer  to  current  transition} 

begin 

for  I:=l  to  NumState  do  begin 
BasicDec[I] ;=0; 

for  J:=0  to  NodeLimit  do  TmpState [J] ;=0; 

CmtDrcupv :  =StateData[I]  ; 
while  (CmtOccupyOnil)  do  begin 

TmpState [CmtOccupy *  .Node]  :=CmtOccupy‘  .Machs ; 

CmtOccupy :  =CmtOccupy *  .  NextOccupy 
end;  {while} 

for  J:=l  to  NumTask  do  begin 
AssnTask [J] : =0 ; 

MaxAssnTask [J] :=0; 
end;  {for} 

CmtTrans  :=TransDataCl]  ; 
while  ( CmtTrans  onil)  do  begin 
J;=CmtTrans*  .TransTask; 
if  (J>0)  then 

MaxAssnTask [J] :=MaxAssnTask[J]+TmpState[CmtTrans* .LossNode]  ; 
CmtTrans ;  =CmtTrans“  .NextTrsuis 
end;  {while} 

for  J:=l  to  NumState  do  TmpDec [ J] : =0.0; 

CmtTrans :  =TransData  [I]  ; 

NextTrans { I , CmtTrans .TmpState .AssnTask .MaxAssnTask , TmpDec .NumTask . 
NumRsrc , NumState . OpsRate .TaskData.RsrcData . NodeData . AltData , Dec . 
BasicDec) 
end;  {for} 

for  I:=0  to  NumState  do  begin 
DecCl.O] :=0.0; 

Dec[I.N\imDec+l]  :=0.0 
end;  {for} 

for  J:=0  to  NumDec  do  Dec [NumState ,J] :=1 .0; 

BasicDec [NumState-t^l]  :=NumDec-)'l ; 

writeln ( Out , ’Number  of  Decision  Variables:  NumDec ;4) 

end;  {MakeDecSpace} 

.{****•********•***•****••***•*•*•**•*•*•**•♦•*•**♦****•*•*•♦•♦•••••**••*•*•••**} 

procedure  OrderAlts(var  NumState:  integer; 

var  AltData:  AltPointer; 
var  Dec:  DecMatrix; 
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var  BasicDec:  StateIntArray) ; 


^ - } 

{  Employs  a  selection  sort  to  order  resource  structure  alternatives  by  } 
{  estimated  performance.  } 
{  > 
{  Calling  Routines:  procedure  Optimize  } 
{  Called  Routines :  none  } 
^ - > 


var  1,1:  integer;  -{local  counter  variables} 

Cmt  Alt,  MinAlt, Last  Alt,  TmpAlt:  AltPointer;  -{pointers  to  local  alternatives} 

begin 

Cmt  Alt :  =AltData ; 

MinAlt :  =Cmt  Alt ; 

while  (CmtAltOnil)  do  begin 

for  J:=l  to  NumDec  do  Dec[NumState-^l,J]  :=0.0: 
with  Cmt  Alt*  do  begin 
CmtEx :  ^Excluded ; 
while  (CmtExOnil)  do  begin 

Dec [NumState-t-i , CmtEx* .ExDec]  :=1 .0; 

CmtEx :  =CmtEx  *  .NextEx 
end;  {while} 

GoodGuess :=0.0; 
for  I;=l  to  HumState  do  begin 
J:=BasicDec[I] ; 

while  (Dec  [NumState-H  ,  J]=l  .0)  do  J:=J-H; 

GoodGuess :  =GoodGuess-»^  1 . 0/Dec  [I ,  J] 
end;  {for} 

if  (GoodGuess>HinAlt* .GoodGuess)  then  MinAlt :=Cmt Alt ; 
end;  {with} 

if  (CmtAlt*  .HextAlt=nil)  then  LastAlt;=CmtAlt ; 

Cmt  Alt :  =CmtAlt*  .HextAlt 
end;  {while} 

while  (MinAltoLastAlt)  do  begin 
if  (AltData=MinAlt)  then  begin 
AltData : =MinAlt  * . Rext Alt ; 

AltData*  .  PrevAlt :  mil 
end  else  begin 

MinAlt* .PrevAlt* . Next Alt :=HinAlt* .HextAlt ; 

MinAlt* .HextAlt* . Prev Alt :=MinAlt* .PrevAlt 
end;  {else} 

TmpAlt : =La8t Alt * . HextAlt ; 

LastAlt* . HextAlt ;=MinAlt; 

MinAlt* .HextAlt :=TmpAlt; 

if  (TmpAltOnil)  then  TmpAlt *. PrevAlt :=MinAlt; 

MinAlt  * . PrevAlt : =La8t Alt ; 

CmtAlt :  =AltData; 

MinAlt :  =Cmt  All ; 

while  (CmtAltoLastAlt* .HextAlt)  do  begin 
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if  (CmtAlt*  .GoodGue8s>MinAlt*  .GoodGuess)  then  MinAlt  :=CmtAlt ; 
CmtAlt :  =CmtAlt* .  NextAIt 
end;  {while> 
end;  {while} 
end;  {OrderAlts} 

■[««*iX«*4t*4i****«*«««**#»**»»»«*»«*»****4i**»***»****»4i»»**i^*»*4i»»«4^*«*4i«[4i*«i^*»***} 

procedure  DualSimplexCvar  NmnState ,NufflDec ,Spec :  integer; 

var  OpMach:  real; 
var  Method:  char; 
var  Opt Alt;  AltPointer; 
var  SpecOpt :  Opt Array ; 

Dec:  DecMatrix; 

BasicDec:  StateIntArray) ; 


^ - } 

{  Determines  optimal  performance  for  a  resource  structure  by  performing  } 
{  dual  pivoting  operations.  } 
{  } 
{  Calling  Routines:  procedure  Optimize  } 
{  Called  Routines:  none  } 
< - } 


var  I,J,Row,Col,NumItn:  integer;  {local  coimter/identif ier  variables} 
Rate, Ratio, Kin, Max:  real;  {local  linear  programming  variables} 

Shortstop:  boolean;  {indicator  for  early  exit} 

begin 

for  I:=l  to  NumState  do  begin 
Col : =Ba8icDec [I] ; 

if  (Dec[NumState+l ,Col]<>0.0)  then 
for  J:=0  to  NumOec-i-1  do 

Dec[HumState+l , J] :=Dec[KumState+l ,J]-Dec[I, J] 
end;  {for} 

Numltn:=0; 

Shortstop : =f alse ; 
repeat 

Row:=0; 

Max; =-0.00001; 
for  I:=l  to  RumState-*-!  do 
if  (0ec[I,0]<Max)  then  begin 
Row:=I; 

Max:=DecCl,0] 
end;  {if} 

if  (RowOO)  then  begin 
Humitn :  =NumItn-t’  1 ; 

Col:=0; 

Min: =99999. 9; 

for  J:  =  l  to  NumDec-t-1  do 

if  (Dec[RoH, J]<0.0)  then  begin 
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Ratio:=Dec[0, J]/Dec[Roi;,J]  ; 

If  (Ratio<Min)  then  begin 
Col:=J; 

Min:=Ratio 
end;  {if} 
end;  {if} 

Rate : =Dec [Ron , Col]  ; 

for  J:=0  to  NnmDec+1  do  Dec [Ron, J] :=Dec [Row, J] /Rate; 
for  I:=0  to  NumState+l  do 

if  (IoRoh)  and  (Dec[I,Col]<>0.0)  then  begin 
Rate:=Dec[l,Col] ; 

for  J:=0  to  NumDec-*-!  do  Dec[I,J]  :=Dec[I,J]-Dec[Row,J]*Rate 
end;  {for} 

BasicDec [Ron] :=Col 
end;  {if} 

OpMach ; =-Dec [0,0]  ; 
case  Method  of 

’N’,’n’:  if  (SpecOpt[Spec]<>nil)  then 

if  (0pMach<=Spec0pt [Spec]* .Optimum)  then  ShortStop:=true; 
’S’, ’s’:  if  (OptAltOnil)  then 

if  (0pMach<=0ptAlt* .Optimum)  then  Shortstop : =true 

end;  {case} 

until  (Ron=0)  or  ShortStop; 

nrite ( Out , ’Expected  Number  of  Operating  Machines:  ’ ,0pMach:8:4) ; 
if  (Ron>0)  then  nriteln(0ut , ’  (or  less)’)  else  nriteln(Out) ; 
nriteln(0ut , ’Iterations :  ’ ,NumItn:4) : 
nriteln(Out) 
end;  {DualSimplex} 

{••••*•••*•**••••***•••♦•••***♦*•**•♦*****•*•***•***••••***•••*•***•**••••*••••} 


procedure  Optimize (var 
var 
var 
var 
var 
var 
var 
var 
var 
var 
var 


NumTask ,NumR8rc,NodeLimit ,NumState ,NumDec :  integer ; 
OpsRate:  real; 

Method:  char; 

TaskData:  TaskArray; 

NodeData:  NodeArray; 

RsrcData:  RsrcArray; 

StateData:  StateArray; 

TransData:  TransArray; 

AltData:  AltPointer; 

Dec;  DecMatrix; 

BasicDec:  State Int Array ) ; 


{ - > 

{  Determines  optimal  resource  structure  using  sequential  linear  programming  } 
{  algorithm.  } 
{  > 
{  Calling  Routines:  program  MAINTOP  } 
{  Called  Routines:  procedure  MakeDecSpace  } 
{ - } 
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var  I, J,RoB,Col,Spec,KumItn;  integer; 

Rate, Ratio, Max, Min:  real; 

TotalCost:  real; 

OpMach:  real; 

Basis:  Stateint Array ; 

SpecOpt:  Opt Array; 

CmtAlt :  AltPointer ; 

Opt Alt:  AltPointer; 

begin 

writeln( ’ **  Optimizing 
Hriteln(Out, >***♦  OPTIMIZATION  SUMMARY  •***’); 
writeln(Out) ; 

NxunDec :  =0 ; 

writelnC  ’  *  Building  Decision  Space _ ’); 

MaReDecSpace ( NumTask , NumRsrc , NodeLimit , NumState , NumDec , OpsRate ,TaBkData , 
NodeData,RsrcData,StateData,AltData,Dec,BasicDec) ; 
if  not  (Method  in  [’1’,’i’])  then  OrderAlts (NumState ,AltData, Dec ,BasicDec) ; 
Hriteln(’  •  Solving  Initial  Linear  Program 
for  Rov:=l  to  NumState  do  begin 
Col:=BasicDec[RoB] ; 

Rate : =Dec [Rob , Col] ; 

for  J:=0  to  NumDec  do  Dec[RoB,J] :=Dec [Rob, J] /Rate; 
for  I:=0  to  NumState  do  if  (IoRob)  and  (Dec[I,Col]<>0.0)  then  begin 
Rate:=Dec[I,Col]  ; 

for  J:=0  to  NumDec  do  Dec[I,J]  :=Dec[I,J]-Dec[RoB,J]*Rate 
end;  -[for} 
end;  -[for} 

Numitn : =NumState ; 
repeat 

Col:=0; 

Max: =0.00001; 

for  J:=l  to  NumDec  do 

if  (Dec[0, J]>Max)  then  begin 
Col:=J; 

Max:=Dec[0,  J] 
end;  -[if} 

if  (ColOO)  then  begin 
Ntunitn :  =NumItn+ 1 ; 

Rob : =0 ; 

Min: =99999. 9; 

for  I:=l  to  NumState  do 

if  (Dec[I ,Col]>0 .0)  then  begin 
Ratio : =Dec [1,0] /Dec [I , Col] ; 
if  (Ratio<Min)  then  begin 
Rob:=I; 

Min:=Ratio 
end;  {if} 
end;  {if} 

Rate : =Dec [Rob , Col]  ; 

for  J:=0  to  NumDec  do  Dec[RoB,J] :=Dec [Row, J] /Rate; 


{local  counter/identifier  variables} 
{local  linear  programming  vsiriables} 
{total  cost  of  resource  structure} 
{number  of  operating  machines} 

{array  of  basic  variable  identifiers} 
{specialty  strategy  optimums} 

{pointer  to  current  alternative} 
{pointer  to  optimal  alternative} 
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for  I:=0  to  NumState+l  do 

if  (lORow)  and  (Dec[l,Col]<>0.0)  then  begin 
Rate : =Dec [I , Col] ; 

for  J;=0  to  NuinDec  do  Dec[I, J]  :=Dec[I,J]-Dec[Row,J]*Rate 
end;  -[for} 

BasicDec[RoH] :=Col 
end;  -tif] 
until  (Col=0); 

writeln(Out Initial  Iterations:  ’  ,)!uinltn:4)  ; 

Hriteln(Out) ; 

writelnC  *  Starting  Sequential  Analysis  _ ’); 

writeln; 

for  I:=l  to  HaxSpec  do  SpecOpt[I] :=nil; 

Cmt  Alt :  =AltData ; 

Opt Alt : =nil ; 

while  (CmtAltOnil)  do  begin 
TotalCo8t:=0.0; 

write(0ut , ’Resource  Structure:  (’); 
for  I:=l  to  NumRsrc  do  begin 

write(Out  .CmtAlt*  .Structured]  :4) ; 

TotalCost :  =TotalCost+CmtAlf  .Structured]  *R8rcDatad]  .Cost ; 
end;  -(for} 
writeln(Out , ’  )’); 

writeln(Out , ’Total  Cost:  ’ .TotalCost :8: 2) ; 
for  J:=0  to  NumDec  do  DecCHuinState-i-i ,  J]  :=0.0; 

DectNumState+l ,NiunDec+l] :=1.0; 
with  Cmt  Alt*  do  begin 
CmtEx :  ^Excluded ; 
while  (CmtExOnil)  do  begin 

Dec  [NumState+l , CmtEx* . ExDec]  :  =1 . 0 ; 

CmtEx :  =CmtEx*  .  NextEx 
end;  -{while} 
end;  -{with} 

Spec :  =Cmt  Alt  *  .  SpecCode ; 

DualSimplex ( NumSt ate , HumDec , Spec , OpMach .Method . Opt Alt , SpecOpt , Dec . 
BasicDec) ; 

CmtAlt*  . Opt innun:=  OpMach; 
if  (OptAltmil)  then  OptAlt  :=CmtAlt 

else  if  (OpMach>OptAlt*  .Optiaw)  then  OptAlt  :=CmtAlt  ; 
if  not  (Method  in  [’S’, ’s’])  then  begin 

if  (SpecOpt  [Spec]  mil)  then  SpecOpt  [Spec]  :  =Cmt Alt  else 

if  (OpMach>SpecOpt [Spec] * . Optimum)  then  SpecOpt [Spec] : =Cmt Alt ; 
end;  -{if} 

writeln( ’Current  Optimum:  ’ .OptAlt* .Optimum: 8:4) ; 

CmtAlt :  =CmtAlt*  .Hext  Alt 
end;  {while} 
writeln(Out) ; 

writeln(Out, ’OPTIMAL  SPECIALIZATIOI  STRATEGY:  ’ .OptAlt* .SpecCode:4) ; 
wr  ite  (Out. ’Resource  Stmcture;  (’); 

for  I:=l  to  NumRsrc  do  wrlte(Out, OptAlt*  .Stmcture [I]  :4) ; 
writeln(Out . ’  )’); 
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vritelnCOut , ’Expected  Number  of  Operating  Machines:  ’ ,0ptAlt~ .Optimum:8 :4) ; 
if  not  (Method  in  [’S’, ’s’])  then  begin 
writeln(Out) ; 
writeln(Out) ; 

writelnC Out, ’OPTIMAL  PERFORMANCE  FOR  EACH  SPECIALIZATION  STRATEGY’); 
Briteln(Out) ; 

I:=l: 

while  (SpecOpt  [I]  Onil)  do  begin 
with  SpecOpt [I]*  do  begin 

writeln( Out , ’Specialization  Strategy:  ’ ,SpecCode :4) ; 
write(0ut , ’Resource  Structure:  (’); 
for  J:=l  to  NumRsrc  do  write( Out , Structure [J] : 4) ; 
writeln(Out , ’  )’); 

writeln(Out , ’Expected  Number  of  Operating  Machines:  ’, 

0ptimiun:8  ;4)  ; 
writeln(Out) 
end;  -[with} 

I:=I+1 

end;  {while} 
end;  {if} 
end;  {Optimize} 


begin 

open(Data, ’MAINTOP.DAT’ ,old) ; 
openCOut,’ MAINTOP. OUT’ ,new) ; 
{assignCData, ’MAINTOP.DAT’ ) ;} 
{assignCOut, ’MAINTOP. OUT’) ;} 
reset (Data) ; 
rewrite (Out) ; 


{VAX/VMS  Pascal} 
{VAX/VMS  Pascal} 
{Turbo-Pascal} 
{Turbo-Pascal} 


writeln( ’RUNNING  PROGRAM  MAINTOP’); 


ReadData ( NumMach , NumTask , NumRsrc ,OpsRate , CostLimit , Method , 
TaskData,RsrcData) ; 
close(Data) ; 

MakeNetwork(NumTask ,NumNode ,0psRate ,TaskData,LowNode ,RodeData) ; 

NodeLimit : =NodeCap(NumMach,NumNode) ; 
if  (NodeLimit<NumNode)  then 

ReduceNetwork(NumTa8k ,NumNode , NodeLimit ,Ta8kData,LowNode ,NodeData) ; 
ListNetwork ( NumTask , NumNode , NodeLimit , NodeData) ; 

MakeStateSpace(NumMach, NodeLimit ,RumState,StateData) ; 

MakeTransSpace (NumTask , NumState , NodeLimit , NodeData , StateData , TransData) ; 
MakeAltSpace(N\imMach, NumTask, NodeLimit .NumRsrc, CostLimit .TaskData, NodeData, 
RsrcData.AltData) ; 

Optimize ( NumTask , NumRsrc .NodeLimit , NumState , NumDec , OpsRate .Method .TaskData , 
NodeData, RsrcData, StateData, TransData, AltData, Dec, BasicDec) ; 
clo8e(0ut) ; 
writeln; 

writeln(’ALL  DONE  (Output  written  to  file  "MAIHTOP.ODT”) ’ ) 
end.  {MAINTOP} 
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B.3  Sample  Output 


INPUT  DATA  SUMMARY 


Number  of  Machines:  2 
Operations  Rate:  0.5000 
Expenditure  Limit:  100.00 
Solution  Method:  NORMAL 

Task  1 

Repair  Rate :  1 . 0000 
Resources  per  Task:  1 
Failure  Rate : 

Prerequisite  Tasks:  2  3 
Task  2 

Repair  Rate :  0 . 2500 
Resources  per  Task:  1 
Failure  Rate:  0.2000 
Prerequisite  Tasks: 

Task  3 

Repair  Rate:  0.5000 
Resources  per  Task:  2 
Failure  Rate:  0.2500 
Prerequisite  Tasks: 


Resoiirce  1 

Unit  Cost:  10.00 

Qualified  Tasks:  1 

Resource  2 

Unit  Cost:  20.00 

Qualified  Tasks:  2 

Resource  3 

Unit  Cost:  25.00 

Qualified  Tasks:  3 

Resource  4 

Unit  Cost:  30.00 

Qualified  Tasks:  2  3 

Resource  5 

Unit  Cost;  33.00 

Qualified  Tasks:  123 


96 


••••  MODEL  SUMMARY 

Number  of  maintenance  nodes  in  full  netsork:  4 


Node  1 

Routing  Probability:  0.5263 

Pending  Tasks :  1 

Eligible  Tasks:  1 

Node  2 

Routing  Probability:  0.1404 

Pending  Tasks :  1  2 

Eligible  Tasks :  2 

Node  3 

Routing  Probability:  0.1880 

Pending  Tasks :  1  3 

Eligible  Tasks:  3 

Node  4 

Routing  Probability:  0.1454 

Pending  Tasks :  123 

Eligible  Tasks:  2  3 

Number  of  system  states:  15 


Specialization 

Strategy 

1 

2 

3 


Resources 
Employed 
{12  3} 
{14} 

{  5  } 
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»»»»  OPTIMIZATION  SUMMARY 

Number  of  Decision  Variables:  35 

Initial  Iterations;  15 

Resource  Structure:  (00003) 

Total  Cost:  99.00 

Expected  Number  of  Operating  Machines:  0.8409 

Iterations ;  3 

Resource  Structure:  (12200) 

Total  Cost;  100.00 

Expected  Number  of  Operating  Machines:  0.8159 

Iterations:  1 

Resource  Structure:  (10030) 

Total  Cost;  100.00 

Expected  Number  of  Operating  Machines:  0.8103 

Iterations :  3 

Resource  Structure;  (21200) 

Total  Cost:  90.00 

Expected  Ntimber  of  Operating  Machines:  0.8080  (or  less) 
Iterations:  3 

Resource  Structure:  (20020) 

Total  Cost:  80.00 

Expected  Number  of  Operating  Machines:  0.8027  (or  less) 
Iterations :  1 


OPTIMAL  SPECIALIZATION  STRATEGY:  3 

Resource  Structure:  (00003) 

Expected  Number  of  Operating  Machines:  0.8409 

OPTIMAL  PERFORMANCE  FOR  EACH  SPECIALIZATION  STRATEGY 

Specialization  Strategy:  1 

Resource  Structure:  (12200) 

Expected  Number  of  Operating  Machines:  0.8159 

Specialization  Strategy:  2 

Resource  Structure:  (  10030) 

Expected  Number  of  Operating  Machines:  0.8103 

Specialization  Strategy:  3 

Resource  Structure:  (00003) 

Expected  Number  of  Operating  Machines:  0.8409 
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